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FOREWORD 


The Industrial Mathematics Society was formed as a forum for 
discussion of matters intermediate betweenthe advanced theory called 
applied mathematics and the routine formulae found in engineering 
handbooks The broad interest in this intermediate field, not only 
among engineers but also among mathematicians has been a source 


of surprise and gratification to the Society. 


In view of this wide interest it seemed desirable to make a per- 
manent record of the papers presented before the various meetings 
of the Society. While there is some hesitancy in starting still another 
scientific journal, it appears that there is need for a publication 


in the intermediate field mentioned above. 


This first issue is a compilation of papers presented during the 
Society’s first year and it is intended to continue on an annual 
basis, at least until the quantity of papers presented before the 


Society forces more frequent publication 
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THE 


The basic principles 





BASIC PRINCIPLES OF CLASSICAL MECHANICS 





of Classical Mechanics may be divided into four categories. 


These categories and the corresponding basic principles are as follows: 


I The 


bran 





Natural Laws, which are universal in character. These laws have 
been developed gradually since man first concerned himself with that 


h of science now known as Mechanics There are four such laws: 


A 


The Law of Conservation of Mass states that mass can neither 
be created nor destroyed Modern physics takes exception to 
this point of view for those processes involving atomic fission, 
however this development is ruled out of Classical Mechanics. 
The Law of Conservation of Mass may be stated mathematically 
as 

d 


—(M = 
at ) 0 


where (M) is the total mass contained in a closed system and 
(t) is time 


The Laws of Motion are associated with the great scientist Sir 
Isaac Newton Modern physics also takes exception to this bas- 
ic principle in the form stated below in that it makes suitable 
modifications. These modifications become important only when 
the velocities of the bodies approach the velocity of light, a 
phenomenon which rarely occurs in Classical Mechanics. The 
Laws of Motion may be stated in vector form as 

i = - 

—(Mv) = F 

dt 

-e 

where (v) is the velocity vector and (F) is the resultant force 
vector acting on the mass (M). 


The Law of Conservation of Energy is frequently called the 
generalized First Law of Thermodynamics. This law, which 
is also invalidated by atomic fission, may be stated mathe- 
matically as 


where (€) is the total energy in a closed system. The total 
energy may be made up of kinetic energy, heat, potential en- 
gy, et 


The Law of Inadmissibility of Decreasing Entropy is also called 
the Secoad Law of Thermodynamics. This law may be stated 
symbolically as 


d 
anil P4 
T°?) 0 


where (S) is the total entropy in a closed system. The use- 
fulness of this law is somewhat limited in those instances in 
which the entropy is not constant. Since this is a common 
situation in nature, this fundamental principle sometimes makes 
it difficult to draw conclusions from our mathematical analyses. 





II Kinematic Relations fall into two broad classes, namely, definitions and 
continuity conditions 


A The definitions needed to formulate the whole of Classical Me- 
chanics in a concise fashion are almost innumerable. A typical 
example of these is the definition of strain 


Ou 

e -=_ 

. Ox 
where (€z) is the elongation per unit length in the (x) direction, 
called strain, and (u) is the displacement in the (x) direction 


B. The continuity conditions insure that the masses being studied 
do not perform in a manner prohibited by the design of the 
System. Examples of this type of conditions are the well-known 
compatibility equations of the theory of elasticit 


ITI Equations of State relate the variables which describe a static state of 
the body to each other They do not involve time Examples of these 
equations are: 


in which (P) is the pressure, (e) is the mass density, (T) is the tem 
perature, and (R) is a constant 


Hooke’s Law for an elastic rod 


where (€) is the strain, (ao) the stress, (T) i the temperature and (E) 
and (a) are the constants 


IV The Dissipation Relations describe the manner in which energy is trans 
formed. These relations involve time and/or time rates of chang [wo 
examples of these relations are: 

Newton’s Law of viscosity for a fluid: 

- iz t ed 'u 

oy 


where ([xy) is the shear stress in the (xy) plans (u) and (v) are the 
velocities in the (x) and (y) directions, respectively, and (mw) is a con 
stant 


Fourier’s Law of heat conduction: 


where (9x) is the rate of heat flow per unit area in the (x) direction, 
(T) is the temperature and (K) is a constant 


It should be noted only the natural laws are of universal application and that the form 


of the remaining basic principles depend upon the problem and material at hand It 
is quite apparent that all of these principles with the exception of equations of state 
involve space and time derivatives For this reason, most of the fundamental e 


quations of mechanics are differential in character 












’ 
A few examples, illustrating the principles stated above and a standard approach for 
continuous bodies are considered below The method utilized, namely that of con- ' 
sidering an infinitesimal element of the body of the form of a rectangular parallele- " 
piped is not the most elegant procedure that can be devised, but it has great sim- 
plicity 
' 
l Example of the Law of Conservation of Mass: 
THE CONTINUITY EQUATION FOR FLUID FLOW 
Consider an infinitesimal volume element in the body as shown. L 


























The Law of Conservation of Mass requires that 

“increase in mass content” = “net inward mass flow” 
Then, if (¢@) is the density and (vx), (vy) and (vz) the velocities in the three co- 
ordinate directions, we have, per unit time: 


3% 


do, 
~dxdydz = dydz - + —— dx) dy 
It xday az ov Gy a2 o(% ax x) y dz 


ov. 
+ pvydzdx - o(y + = Y ) da dx 
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Ov, ; 
+ ov,dxdy - o(vy, + ——dz) dxdy 
as 
or 

dp . oy av av 

= @& @—9( a= ¢ at + a 

dt Ox y 02 

where ($f) is the total rate of change of the density. This is the usual form 
t 


of the continuity equation 
2 Example of the Law of Conservation of Energy: 
THE HEAT CONDUCTION EQUATION 
Consider an infinitesimal volume element in the body as shown, 
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j The Law of Conservation of Energy requires that 
“increase in energy content” = “net inward energy flow” 
; Then, if (c) is the specific heat, (?) is the mass density, (T) is the temperature 
and (9, ), (4%) and (gz) are the rates of heat flow in the three coordinate di- 


rections, we have, per unit time: 


oT , ‘ oq , 
Co—dxdydz = q dydz - (q. + —&dx) dydz 
Ot . - Ox 


oq 
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— og finn 
+ q.dxdy - (q, + = dz )dx dj 
or 
oT 
Co— S «= ( 04. + oq, + O49, ) 
ot Ox Oy Oz 


It is common to introduce Fourier’s Law of heat conduction into this equation 
and thus obtain the usual form: 





aT zs, .38. 27 
ot Co Ox oy 02 


3 Example of the Laws of Motion: 
EULER'S EQUATIONS FOR FLUID FLOW 
Consider an infinitesimal volume of the fluid as shown 
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The Laws of Motion require that 


“rate of change of momentum” = “resultant force 


4 oP . . 
== (ov, dxdydz) = Pdydz -(P+ =~ dx )dy dz + oX dx dy dz 


where (pdx dy 42) )== the mass of the element, (P) is the pressure and (X) is the 
body force in the (x) direction per unit mass. Thus, for each direction, there 
is one equation of the form 





dv,  ~ r 
eit ~ 3x fn 
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where (<¥) is the total time rate of change of the velocity. 
































+ Example of the Laws of Motion 
aTcormnrr ™ , an ro * , ~ + . ; 
LONGITUDINAL VIBRATIONS OF AN ELASTIC BAR 
Consider an infinitesimal slice of an elastic bar as shown. ; 
eS ee ee aS a ae 
z “ 
“—-—-=-- “|——-—~+1 
| . | Ta ag, tPF 
Oo + <2 dx 
x | 
; 
in—_ Gus cas Giee deme oes op a ee a ae 
a 
| 
@) ———___—— x ; 
The Laws of Motion require that 
“rate of change of momentum” = ‘resultant force” 
3G, Oa 
——-(pAVdx)= -gA +( got—adax)A 
ot Ox 
where (p) and (A) are nearly constant for small motions, and (oc) is the tensile 
stres in the bar Thus 
OV Oc 
= * = 
ot Ox 
It is common to introduce the definition of velocity and strain as well as Hooke’s 
Law 
Ou ou - 
ys -- ¢ = os oc be 
ot Ox 
where (u) the displacement in the (x) direction. The basic equation may now 
be written as 
» 8 ' 
ou E 97u 
ot? re) ox? 
which the usual form of this equation 
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How can we protect against this calamity? Promote if you must, but plan on the 
bas of detailed studie not on promotior Let the experienced construction man 

an the program, Force the con te ientist to develop all deta before design 
and construction, r the latter will be a true gambl« 
During the progress of the work beware of new proposa Changes toward improve 
ment rr more dangerously, to ex efulne of the computer generally set 
back the overall program. It is diffi t to avoid being persuaded to add a shifting 
register which would only add one panel of 50 tubes and can be placed in the space 
at the top of the input cabinet so that aerodynamic problems for ultrasonic speeds 
can be solved in one-third of the time otherwise required For example, our EDVAC 
was originally planned to have 200 tubes In 1947 it was up to about 4, 000 tubes, 
Through change to crystal diode gating rcuit t dropped to 2,000 tubes and 5,000 
crystal diodes. By early 1948 it was almost 4,000 tubes and 8,000 crystal diodes. 
True, the final EDVAC was far better than that planned with 1, 200 tube but a less 
satisfactory computer available for problems much earlier might well have been more 
valuable 
Let us build some computers but refuse to expect a machine under nstruction t« 
keep abreast with the currently new proposals in the field 
Reliability of Digital Computer 
In the EDVAC, information, be it number or instruction, 1 1 SuCcCE ion of go or no 
go conditions, the presence or absence of pulses, In many of these rcuits pulses, 
or the absence of pulses, are passed once every microsecond, Should a tube or a 
crystal rectifier miss a beat, an error almost certain to occur Thus, with the 
number of such elements in the EDVA( ne missed pulse in 1012 would cause an 
error about once an hour, A probability of error of one ten billionth of one per cent 
would mean that problems of an hour duration would be in trouble. 
To transmit audible sounds of extremely high quality requires about 100,000 yeas or 
nayes per second, An error of one of these would correspond to the lk of one-half 
cycle of the 10,000 cycle component pe: econd about the equivalent of one grain 
of dust on a very high quality phonograph record 
Thus the reliability of computer I t be some 10 n n meé is good a in 
excellent quality audio amplifier. Presented in th Nay mputer reliability ap 
pears to be a staggering problem. Howeve1 he tu n in fact quite le evere 
Let us consider first the random err« hose error where Ise ost « 
gained by chance due t: truly unordered random current ich as that of the pitter 
patter of electrons triking the plate of lum tube Here the obal ty of an 
interfering pulse being irge enougn t ympare with the niormatior rie i an 
exponential function of the amplitude Because of tl exponent haracte of the 
probability, we in Say loosely to be ‘ that it twice as | ito ge yroba 
bility of error f one in 10/4 t get ne u o* Ever t the 
present time, the probability of herr ¢ ha ne in 10/4 
Were we todo half again as good - exponential we V have a probable 
error of less than one in one hundred ye 

A more areful analysis of this type f error indicate that I expe tion 
very reasonable However, relaxing our effort it gh ma AuS¢E n ibearable 
number of errors t« e encounters In many ice s the n nen proacl 
failure, the random error ncrease e r i Some tem « reventative 
maintenance is high essentia 
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Error jue to failure of the rcuitry are aln alway I 
rdered ert Here eve! tne ry lest necking yvstem 
re iit | ive 1id I t alWay | iG ses are tn 
ystemat it hidden error F or nstance ertain 
tunct yn niy m tne ise f the east lgniti in lgit: i 

tent but hidden error I believe that much in e dor 
the I iter t avoid ich |} ties, (Thi nsider 
idvantage f the er i type mas ne n ¢ rast t tne f 
The need for automat hecking ystems ha been th S 
However these have eer irge d sions and tudl 
Want t ra e tne juest nm I the id l if t I any 
Checking ysten vt h are intended t find random e 
effect hould be made unnecessary by somewhat more 
vith a preventative maintenance progran Simp »CcKI 
roblen needed t ‘ part failure fa 
hould be attacked by considered design, 
Preventative maintenance directed toward keeping the jf 
ifely be A tne danger point 11Ti 1it Tt i omp ish 
ilacement t il] mponent ict is tubpe hich rit 
would be iff ient 4 ntere ng roy ul te 
r other parameter f the design whict oO “¢ 
ranaom error ource 
Checking ysten generally require nsiderable equi 
cause of trouble Then again hecking systen isef 
t finds very many ch error the usefu Omputing tin 
he< ng ysten ire iluable if they find a few errors 
iseful if the checking system were t find many error 
immary 
Digital mputer ire wonderf machine part i] A 
potentiality, There r lerable danger of over-pron 
i rogran hould be done n a detailed basis and not on 
pronouncement 
Reliability 1 ’ roblen f an rder new to electronic 
ment t late ind ate that the required li lity is no 
Recent cont ms that elaborate hecking sy 
exarr nea ince they ire I m ise lI the. 
whole a fa re if they find many error 
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“ANALOGIES IN FLUID FLOW, ELASTICITY AND SOAP FILMS” 
Dr. G. E. Hay! 


l Introduction. In this paper it is proposed to present some of the fundamentals 
of the mathematical theories of elasticity, hydrodynamics, and soap films. Each of 
these three theories is quite extensive, and can easily constitute the subject matter 
of a large number of papers such as the present one. Accordingly only enough of 


the fundamentals of these three theories will be presented here in order to demon- 
Strate an analogy between them. As we shall see, this analogy exists because three 
problems, one from each of the three fields mentioned above, all reduce to the same 
problem in mathematics These three problems are the following: (a) the torsion 
of a beam of uniform cross-section; (b) the two-dimensional flow of an incompressible 
perfect fluid inside a hollow rotating cylinder; (c) the equilibrium of a soap film 
stretched across the open end of a tube, the pressure inside the tube being greater 
than that outside the tube 


The importance of the analogy under consideration here lies in the fact that it is 
possible to obtain information about the solutions of any two of the above three prob- 
lems through experimental work carried out on the third problem. 


ELASTICITY 





2. <A Study of Stress A body is said to be elastic if it suffers no permanent dis- 
tortion when a load is applied to it and is then removed. This property of being 
elastic then depends both on the nature of the body and on the magnitude of the load. 
In this paper we shall consider only elastic bodies 





Let us consider a general elastic body, such as the one occupying the region (R) in 
Figure 1 We shall neglect effects due to gravity. In order to describe the behavior 
of this body mathematically, we introduce rectangular cartesian coordinates (x)), 
(x9), (x3), as shown Throughout this talk we shall adopt the convention that Latin 
subs ripts range over the values 1, 2, 3 Accordingly we can denote the three 


rectangular cartesian coordinates collectively by the single symbol (x;). 


eK, 
NAS 
R J 
— tee lea tle Xs 
ge Figure 1 
We shall now introduce functions which describe the internal reactions between all 
pairs of adjacent particles in the body Let (P) be a general particle in the body, 
with coordinates (x;). Let (1;) be the components of a unit vector at (P), and let 
s 


(A) be the plane through (P) and perpendicular to the unit vector (lj). The stress 





Dr. G. E. Hay is in the Mathematics Department of the University of Michigan. 
This paper was presented at a meeting of the Industrial Mathematics Society on 
Friday, November 4, 1949 











, “ ; 
vector (T;) can now be defined It is the force per unit area at (P) exerted across 
i i 
(A), by those part é which lie on the same side of (A) as the vector (1;), on those ; 
partic le on the other le of (4) Obviously (T;) are three functions of the co- 
ordinate (x of (P) and the omponents of the unit vector (l;) normal to (4). To 
lienote tr At an write 
' J 
(2 I | (x 
We ha now introduce the tre omponents There are nine of them, which we 
lenot t the ymbol (T The three omponents (Tj ;) are the omponents of the 
tre vector at (P) when (1l;) made to point in the direction of the positive (x,) 
ixi The definitior f (T5;) and (T..) are similar 
é) 3} 
, 
“ Figure 2 
Let us now cor ler the tetrahedron (T) shown in Figure 2 One corner of (T) is at 
the oint (P), and the three fa through (P) are parallel to the oordinate planes 
The fourth face paralle to the plane (4) The total force acting across the four 
face must be qual to zer Since the elastic body is in equilibrium This total 
force an be expr‘ 1 in tern if the values of (T;) and (T;;) at (P), as well as 
(l;s) and the area of the loping fa orrect to the lowest order in the dimensions 
of (T) If (T) then made to llapse on the point (P), the following relation re- 
ilt 
(2. 2) , 227 
In order to determine the reactions between all pairs of adjacent particles in an e- 
astic body, it is thus necessary only to find the nine stress components (T;;), which 
. J 
are functions of the oordinates 
Let us now consider the box (B) shown in Figure 3 One corner of (B) is at the 
point (P), and all faces of the box are parallel to the coordinate planes Since the 
elastic body is in equilibrium, the forces acting across the faces of (B) must have 
a zero sum anda zero moment about the origin (O) of the coordinate system. These 
forces can be expressed in terms of the values of (T..) at the point (P), orrect to 
; . 1) : ad 
the lowest order in the dimensions of (B) If (B) is° then made to collapse on the 
point (P), the following equations result: 
3 
(2.3) ~~ at 
= ee w O, 
j=l Oxj 
(2.4) I I 
The three equations in (2.3) are called the equations of equilibrium We note that 
because of (2.4) it is mecessary to find only six stress cemponents in order to de- 
termine the reactions between all pairs of adjacent particles in the body 
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Figure 3 
3 A Study of Strain Let (uj) denote the three components of the displacement of 








the general particle (P) in the elastic body. The strain components (ejj) are then 
defined by the relations 





au du 
ee | 
¢ 1 > - 1/ + : 
(3.1) eij = 1/2 ( 3x; oxi ) 
buc Sirail COMponents Mave physical Meanings, ror exampie, (|e ) 1S tne increase 


in length per unit length of an element at (P) originally parallel to the (x,) axis; 
also (C53) is one half of the decrease in the angle between the elements at (P) o- 
riginally parallel to the (x9) and (x,) axes 


It is interesting to note that (e;; O) identically when (u;) are the components of a 
rigid-body displacement 


In many problems it is convenient to find (€;j) first, then obtain (uj) from equations 
(3.1) These equations (3.1) are then regarded as six equations for the determi- 
nation of the three unknowns (u;) These six equations possess solutions only when 
(ejj) satisfy a set of six second order differential equations called the equations of 
compatibility. These equations will not be required in our discussions which follow, 
Accordingly they will not be reproduced here. They are mentioned only for com- 
pleteness 





4. The Stress-Strain Relations. Since stress and strain ge hand-in-hand, there 
must be some relation between them. The usual form of this relation is given by 
a law called the Generalized Hooke’s Law which is as follows: the stress com- 
ponents are linear functions of the strain components. 











We shall assume that the elastic body under consideration is isotropic; this means 
that all directions in the body are equivalent elastically. In this case the stress- 
Strain relations take the form 


(4.1) = nO5 i; + 2 ej 


1) J 


where ()) and (xz) are elastic constants, and 








3 
(4. 2) O = 2. kk 
k=1 ° 
(4.3) d ij =! if i = j, 
/ 
= O 2s Sa 
9. The Torsion of a Circular Shaft Let us consider a shaft with a uniform circu- 
lar cross-section, the shaft being acted upon by equal and opposite twisting couples 
(M), as shown in Figure 4. In order to determine the stress and strain in the beam 


bility for the unknowns (Tjj), (€;;), (uj), where the solutions must be such that the 
proper loading conditions are satisfied over the sides and ends of the shaft. 


t is necessary to solve equations (2.3), (2.4), (3.1) and the equations of compati- 








| 1 ( Xa 


\ straight forward attack leads to difficulties, because of the large number of un- 
knowns and the complexities of the equations involved Accordingly it is customary 
to mak rtain simplifying assumptions, then solve the problem, justifying the as- 
imptior by an appeal to the fact it is known that the mathematical problem has 
but ne olution rhe implifying assumptions made in the present case are the 
following a) particle are not displaced out of the cross-section in which they lay 
nitially, (b) each cro ection i rotated as a rigid body through an angle pro- 
portional to the listance from one end of the shaft On the basis of these two as- 
imptior t the readily follow that 


it I juat ju b r se we the aiding ynditions over 
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f The Tor no Shaft ral Cross-Section Let us now consider the 

prob em of t rsio 1 5, but in the case when thi ross-section is 

1 general region (S) bounded by a curve (C), as shown in Figure 5. In order to 

find the tre and strain in the shaft, implifying assumptions are made just as in 

Section 5 In the present ase these assumptions are: (a) the omponents of dis- 
acement (u;) and (us) are given by (5.1); (b) the axial omponent of displacement 
la) a funct f (x,) and (x5) only We an thu te in the resent cas¢ 


Because of (6.1 the trail mponents (¢ ) can then be computed in terms of (@) 
from (3 From (4 the stre omportents (Tj ;) can then be found in terms of 
~ ) These tre omponent atisfy the equations of equilibrium (2.3) if 
) ) 
e 9 a<¢ a“? 
{ Zz) > 
x" . 3x Q in » 
) 
These tre ompor i ( it y t oading ondit ns on t 5 5 I tr snatt 
roy ] ] 
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Figure 5 


where (mj;) are the components of the unit vector normal to the sides of the shaft, 
and (d@/dm) denotes the directional derivative of (@) in the direction of the vector 
(mj). Just as in the case of the circular shaft, these stress components satisfy the 
loading conditions on the ends of the shaft only when the twisting couples (M) arse 
distributed over the ends of the shaft in a special way We then have the following 
onclusion: in order to determine the stress and strain in a beam with a general 








ross-section and in a state of torsion, it is only necessary to find a function (P) 





satisfying the differential equation (6.2) throughout the cross-section (S) of the beam, 














and the boundary condition (6.2) on the boundary (C) of the cross-section, Of course, 
in general this solution will be valid everywhere except near the ends of the beam, 
Any function satisfying the differential equation (6.2) is said to be harmoni« Every 
harmonic function possesses a conjugate function which is also harmonic and is re- 
lated to it by certain relations called the cauchy-riemann relations If (>) denotes 
the conjugate of the warping function (#), these relations ar« 
(f 4) 3 3 4 aw 
aX] xD ax] 
Because of (6.3) and (6.4), we see that the conditions for the determination of (1 ) 
ire now 
32% 4 27 
(6.5) as a O in $ 
ax* ax> 
2 
9 9 
(6. 6) y 2 (x, + X5 ) + const on C 
The solving of the torsion problem can thus also be reduced to the determination of 
a function (¥) satisfying the differential equation (6.5) throughout the cross-section, 
and boundary condition (6.6) on the boundary of the cross-section This is so be- 


ause the warping function (@) can be readily deduced from (6.4), once (y) has 


been determined 


There is a third mathematical formulation of the torsion problem which is of interest 
We introduce a function (y) by the relation 
2 2 
(6. 7) X= w-1/2(x, + x,) 
| - 


It is then readily seen that the solving of the torsion problem can be reduced to the 


etermination of a function (X) satisfying the equations 


9 2 
(6. 8) ot al in »S, 
ox + ax, 2 


(f 4) V const on C 








HYDRODYNAMICS 








7 Irrotational Flow of a Fluid Let u ynsider a mass of incompressible fluid in 
motion Let us introduce rectangular cartesian coordinates (x;), as before. Let (P) 
be a general particle in the fluid, and let (v;) denote the components of the velocity 
vector of (P) In general (v;) are functions of the coordinates of (P) and of time 1 
(t) [To denote this we writs 
(7.1) Vv; ; (x, t) 
This flow is said to be irrotational if the following relations are satisfied identical] 
(7. 2) 2 2) O 
o*j o*i 
The study of irrotational flows i isually simpler than the study of rotational flows, 
ince inthe former case the determination of the three functions (v;) can be reduced 
to the determination of a single function called a velocity potential 
Two-Dimensional Flow A flow i said to be two-dimensional if it is possible 
to choose coordinate axe o that 
a" y*2 
{ = V ) 
Ax 0X2 
o ; } 
Let us now consider a two-dimensional flow of an incompressible fluid Further let 
the flow be irrotational; th means that 
v9 
(8 2) 
AX9 »x 
In addition, let us suppose the flow to be steady; this means that 
(¢ 3) Vi Vi (Xy Xs ) V9 = V5 (X} X95 ), 


ie.,the distribution of velocities throughout the fluid does not change with the time 


9 





Let (P) be a point in the (x;X»5) plane and let (C*) be any curve joining (P) to the 
origin (O) as shown in Figure ¢f Let us consider the cylinder made up of all the 
straight lines through (C!) and parallel to the (Xq) axis The amount of fluid flowing 


across a unit length of this cylinder per unit time, in the direction of the arrow, 


is denoted by (‘\) Since the fluid is incompressible (W) is a function only of the 
coordinates (x;) and (x5) of (P) It is.called the stream function The lines (1) = 
const. are called the stream line 





now change the corrdinate (x,) of (P) by an amount (dx,;), as shown in Figure 


f This produces a change of (d¥y) in (\), which change is due to the rate of 


flow 
across a unit length of the cylinder with (PQ) 


as its trace This rate of flow is 
easily seen to be (v9dx,) whence we have 


(xg) of (P) to change by an amount (dxg) we find 


(8. 5) v, © - FY 
0X2 


[The determination of the two velocity components has thus been reduced to the 


de 
termination of the single stream function (W) 


Because of (8.2), (8.4), (8.5) it follows that 


9 9 
- = , 29 O 
(8.6) 
9x; > x. 
l ) 
9 Motion of a Fluid Inside a Rotating Cylinder Let us now consider the steady 
irrotational flow of a perfect fluid inside a hollow cylinder which is rotating with a 


iniform angular velocity («) about a line parallel to the generators of the cylinder 
Let (S) denote the cross section of this cylinder, and let (C) denote the curve bound 


A ‘ 


ing the cross section, as shown in Figure 


The motion is characterized by a stream function (|) which satisfies the differential 
equation (8 6) throughout the region (S). It remains only to deduce the boundary con 


dition (W) must satisfy on the boundary (C) of the cross section 


Let (P) be a general point on the boundary (C), as shown. 


The velocity of (P) is 
(Wr) in the direction shown, 


and the component of this velocity in the direction normal 
to (C) is (wr cosa), where (&) is the angle between the radius vector to (P) and the 
tangent to (C) at (P), as shown. Now the velocity in the direction of the normal to 
(C) of the fluid particle in contact with (P) must also be (Wr cos&). Hence the rate 
of flow across a unit length of the cylinder with the element (PQ) of (C) as its trace 
is (Wr coseds), where (ds) is the length of (PQ). But this is equal to (d}). Thus 


(9 1) av 


Wr cOSe& on C 
2 c 
Bun 
From Figure 7 it fe lows that ( o or) Thus 
ds 
a i . 
¥ = wr dr on € 
A s ds 


whence we have the boundary condition 





; 2 ; 
(9. 2) W = 1/2 wr” + const on ¢ 
= y 
, : yA 2 
Now r“ = 1/2 (x; + X9) If we introduce new rectangular cartesian coordinates (Y¥}), 


(Yo) defined by the relations 
wa 


Y) =NW xy 3 > =¥w xo 








We then see that the problem under nsideration is characterized by a function (YW) 


which satisfies the differential equatior 
) , 
» 1 
| E>, } ¥ ¥ () nn» 
Aye AUZ 
} 
} P 
) 4 a’ i=. * Sh 4 nst n ¢ 
) 
It hould be noted that equations (6.5) and (6.6) for the determination of the torsion 
inction (WW), and equations (9.3) and (9.4) for the determination of the stream function 
(WY), are mathematically identical Thi jemonstrates one analogy between elasticity 


10 The Equilibrium of a Soap Film Let us consider a thin-walled tube, one end 
yf which is closed by a ap film Let the region occupied by the soap film be the 


region (S) shown in Figure 8, the boundary of (S) being a curve (C) Let us now 
apply a pressure to the air inside the tube, denoting by (p) the excess of the pressure 
nside the tube over the pre ire itside the tube Of course the soap film bulges 

itward We shall now onsider the equation for the determination of the shape of 
th soap filn 


p 
. Figure 8 

Let (T) denote the surface tension of the soap film It is the force per unit length 

acting across any curve drawn on the soap film, and is a constant for any given 

oap tilm 


Let (P) be a general point in the undisplaced soap film, and let the coordinates of 
(P) be (x;), (X9) Let (W) be the displacement of (P) normal to the region (S) It 
then required to find (W) as a function of (x,) and (x9) 


Let u onsider an element of the undisplaced film at (P), with edges of lengths (dx,), 
(dx9), a shown This element is displaced into a curved surface in equilibrium under 
the pressure (p) and irface tensions acting along its four edges The force acting 
on the element in the direction of the (X3) axis and due to the pressure (p) is (pdx ,dx9) 
The resultant in thi same direction Af the surface tension forces a ting on the two 
edges originally parallel to the (x,) ajc can be shown to be 


I ; GX)aX5 
>. 2 y, 
dx, 


when the assumption i made that the displacement (W) 1s much smaller than the 


4 


smallest diameter of (S) Similarly the resultant of the surface tension forces acting 


n the two edges originally parallel to the (x5) axis can be shown to be 


2° W 
I = ix, dx9 
Ox - 





that the total force acting on this element in the direction of the 


o W i 

+ 7 ty —— 

> x “ Ox‘ I 
Ux, 


equilibrium of a soap film under small deflection 


e new rectangular cartesian coordinates (2 1)» (Z.) defined by the 


< 


r 1 
= x Zz > 4 X49 
V2 ” ¥ 27 ” 
1 (10.2) for the determination of the deflection (W) of the 
€ IT i 
) >” W . 
, : 2 n §, 
12 as © 
i a 
() on ¢ 
ne equations (6. 8), (6.9) for the determination of the torsion function 
+), (10.5) for the determination of the deflection (W) of a soap 


unt in equation (6.9) can be chosen equal to zero with no loss 


appears that the problem of determining the torsion function and 


ithematically identical This demonstrates an analogy between 
ve analogy, the following statements can be verified: 
ne swept out by the inflated membrane is proportional to the torsional 


responding shaft in torsion, where the torsional rigidity is the 
equired to produce a unit angle of twist per unit length Because 
imple determination of the comparative torsional rigidities of a 


be made by experimental work carried out with soap films. 


the displaced soap film where the gradient is steepest corre 
visted shaft where the stress is greatest 


ve been constructed for measuring and interpreting the per- 
nm n order to thereby obtain information about the correspond- 


e ticity and hydrodynami 


hould be pointed out that the analogies discussed here are not the 
between these three fields Further information concerning these 

1ined from the following books: H. Lamb, Hydrodynamics, Cam- 
re Cambridge, England, 6th ed., 1932; S. P. ~ Timoshenko, 


McGraw-Hill Book Co., New York, 1934 
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“On the Importance of Anchorage in the 
mae 


Applications of Mathematics to Industry 
l 


Richard S_ Burington 


Introduction There is a growing bond between the methods used by men in 
management and those used by scientists. The potential possibilities of this bond are 
o great that the power of such developments should be more fully understood. This 
is a matter of concern for executives, educators, scientists, and engineers. In par- 


ticular, those engineers, scientists and mathematicians whose interests lie in border 
fields or in industry should be especially concerned 

Most of us are well acquainted with the contributions scientists, engineers and mathe- 
maticians have made in the discovery, research, development, production, and the 
realization of all sorts of industrial achievements; the improvement of transportation-- 
by land, by sea, and by air--, the wide-spread introduction of labor-saving devices, 


and the great material progress of recent years 


Many spectacular achievements have been widely publicized; while relatively little is 


known by the public at large concerning the less spectacular contributions professional 
men ich as applied mathematicians have made These contributions have been 
wide spread Each of us i familiar with some specific phases of this progress 

Perhaps less is known by the publi and even by the professions themselves, of the 
contributions professional men--and, in particular, applied and industrial mathema- 
tician can ard have made to the planning, execution, coordination, and administration 
of large-scale enterprise, and developments and operations of every conceivable type. 


These contributions, whether they have involved purely technical research, develop- 


ment, and analysis in the small, or whether they have entailed analysis relating to 
the planning, execution and administration--that is management--of overall large-scale 
developments or operation have in every successful case hinged upon studies which 


were firmly anchored in the fields involved 


As with the problems of nature, the problems of the industrial mathematician are 
eldom staged in the manner mathematicians might wish they were The problems 
the industrial mathematician must face are commonly those in which a portion of the 


premises and perhaps a portion of the conclusions are given; and to varying extent 
portions of the premises and conclusions may be given as data with different degrees 
of probability rather than as precisely determined data To form well - mapped 
structures joining the assumptions and conclusions, rigorous logical reasoning is 
needed, but it is the foundation upon which the structure rests, the manner in which 
it is anchored, which is most important, and which is so often uncertain and hence 


precarious 


The great importance which must be attached to the proper anchorage of analysis to 
the underlying fields of investigation cannot be over-emphasized This is particularly 
true in the applications of mathematical and statistical methods to the problems of 
industry, both in the technical as well as in the management fields 





1 Dr. Burington is Director of the Evaluation and Analyses Group and Chief Mathe- 
matician, Bureau of Ordinance, Navy Department. The views expressed in this paper 
are those of the author and are not necessarily those of the Navy Department This 


paper was delivered at a meeting of the Society on May 24, 1949, 
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2 The formulation of a problem, a fundamental process. Most of you have had the 
experience of being confronted with stidies of an industrial, physical or natural 
character in which the confusion was almost overwhelming, and have found it extremely 
difficult to extract from the mass, or dirth, of available information a sensible formu- 
lation of the pertinent problems. The larger and more difficult problems of industry 
are apt to be of this type. Just how to determine what is essential and fundamental 
is very often distressingly difficult An overall point of view is usually necessary 
Sometimes dealing with the simplest cases helps. But in all cases a thorough under- 
standing of the underlying physical and economic fields is extremely valuable. Fre 
quently, many attempts are necessary Occasionally, unconventional treatments are 
fruitful. The product of such considerations is a certain simplification and idealization 
of the original situation and formulation of the problem at hand. It is such idealized 
models which are made the subject of analysis 

The use of approximate equivalent systems appears as part of a fundamental process 
in almost all studies of physical, economic and industrial phenomena One method 
of approach consists broadly of the following steps: 


(1) The extraction from the phenomena or situation S of a nearly equivalent 
physical (economic, or other) model P 


(2) Reduction of the physical (or other) model P to an approximate equivalent 
mathematical model M of amenable treatment This may involve the selection 
of a set of observable quantities which are assumed to be appropriate for the 
studies at hand 


(3) A study of mathematical model M involving various relationships R between 
the observable and perhaps non-observable quantities, with solutions of pertinent 
items relating to the system M 


(4) The interpretations of the solutions found in (3) in terms of the mathematical 
model M 


(5) The interpretation of the solutions found in (4) in terms of the physical 
(or other) model P 


(6) The interpretation of the results (5) in the original physical (or other) 
settings S 


Steps (1) or (2) constitute the formulation of the foundation for the mathemati tudie 


involved. If these foundations are not firmly anchored inthe underlying basic setting 
S, then the results obtained in (5) and (6) may be properly questioned 


Steps (3) and (4) are properly those for the mathematician as such. It is his business 


o reveal all that is implied by the relationships R This step may not always be 
possible with the present state of the art. However, the mathematician may then be 
forced to make further modifications of the model M with the attended dangers to the 
soundness of his formulation and anchorage in the original phenomena §S. Such situ 
ations sometimes lead to important researches in pure mathematic 


It is of considerable importance to realize that idealized mathematical model M i 
somewhat different from the actual prototype situation §S For this reason mathe 


matical results derived from the model M and subsequently interpreted inS may and 
Should be questioned as to their validity, applicability and relevance to the original 
setting S. Thus the anchorage of M in S must be the subject of thorough investi 
gation, in the light of its bearing on the relevance and applicability to the original and 





primary setting S of the mathematical results derived from the model. All of us who 
deal with applied and in particular indusfrial and business problems, know that if the 
model differs too greatly from the original situation S, any related practical con- 
iderations and projects may be delayed or brought to naught. No matter how logical 
or beautiful an analysis is made in connection with an industrial consideration, if it 
fails to meet the realistic tests of relevance, adequacy and applicability, sooner or 
later it will be recognized by those engaged in the original larger problems as being 


inadequate and impractical and may even be called worthless 


[To determine whether a mode! falls within allowable bounds must, as a general rule, 
be determined by experiment, after the analysis or at least a portion of it has been 
considered. One crucial test as to the worth of a theory or design, deduced mathe- 
matically or otherwise, hinges on the degree of agreement or disagreement at suitable 


points with experiment und measurements This usually entails a consideration of 
the reliability of the test data and the reliability of the premises upon which the 
inalyse were made Thus asSuming the absence of actual logical or calculating 
errors, and the legitimacy of the pure mathematical processes used, the inferences 


to be drawn from the final judgment of the worth of a theory or design bear upon 
the legitimacy of the abstractions and simplifications which were made in constructing 


the model and anchoring it in the original setting S 
Another crucial test of a theory design yr particular model structure lies in the 
ower of the theory to predict observations in the prototype field previously unknown, 
and in the ability f the theory to include all known pertinent facts in the universe 
t lesigned t ie ribe 
If, in addition i particular theory or model structure proves to be suggestive and 
timulating to new development it is very likely to turn out to be a truly successful 
theory and perhay more than some other apparently suitable working theories 
which are le iggestive 
In examining the typical equations developed or encountered in mathematical physics-- 
and in industry ne soon realizes that they relate partly to physical objects and 
event und partly to a mathematical structure which was introduced for reference 
irpose Further study of ich equations will reveal a distribution between those 
juantities called observable which can be determined by a specified observational 
rocedure and those quantitie which are not observables The non-observables may 
depend at least partly on a mathematical framework which cannot be the subject of 
actual observation The non-observables may be very useful for computational and 
theoretical or descriptive purpose However, such non-observables must be removed 


in the final stage where actual predictions must be checked by observational verifiable 


result A consciousne of this aspect of anchorage is most important 


[The importance of experiment as an aid inthe formulation of a well-anchored mathe- 


matical framework for a field of study should not be overlooked. Consider for example 
the old problem of attempting to understand and predict cavitation, say in the neigh- 
borhood of a propeller, turbine blade, or other body immersed in a fluid Many 
attempts t inderstand the growth and collapse of cavitation bubbles have been made 
because of their bearing on the performance of various machines, et< Were one to 
make a new start in formulating ich a theory of cavitation, one should study every 
yit of experimental data available: for example, the reports and photographs made 


it a very high speed of the formation and collapse of bubbles from tests recently 
made in tunnels at the California Institute of Technology From this one should try 
to draw a certain mechanistic picture of the phenomena which could be colored with 
known thermodynamic and hydrodynamic knowledge In this way the construction of 


1 well-anchored structure for an improved theory of cavitation would begin 








It 1S very important to 


experimental science that theories be proposed. Without 
theory much experimentation would be fruitless, or poorly directed As a tentative 
theory is formulated, particular attention to its agreement with observations should 
be made. A slight inconsistency of theory with observations may lead to new dis 
coveries which (without a theory) might otherwise remain unnoticed These dis 


crepancies may lead to an improved theory, 


and so on; but always the anchorage of 


theory to observations must be kept in mind 

The industrial mathematician, scientist and engineer inthe course of his work some 
times becomes concerned with a mathematical structure of a theory for an applied 
field, which is considered in the light of its predictions as a basis for making im 
portant designs, and which may involve considerable financial, personnel or material 
risk In such cases the experienced qualitative or quantitative judgments of other 
should be sought, a judgment which will take into consideration factors quite apart 
from the particular mathematical structures used in the work, a judgment which may 
take account of the reliability of data, the relative importance of various factors and 
other pertinent items beyond the scope of the mathematical structure available, each 
with appropriate weights In making final decisions--as must be done in industrial 
problems-- areful synthesis of all of the approaches and method available should 
be made to ensure the best possible decision 

3 Examples strating the importance of trong anchorage I have ketched 
an outline of a fundamental proce fan ar to all of u the proce f formulating 
a working theor r model for the r e of studying some field, physical or other 
wise; and in industry often for the rpose of developing producing or improving 
1 product, a proce in operation yr the ke Thi art of formulating a problem 
na specif field of application isually the more difficult and often it more 
mportant than it actual solution 

In order t« trate this proces nd to di ISS in more detail the importance of 
good anchorage in the study of industrial of other problem I shall di J briefly 
ertain example which are suggestive of the importance of the theme of this lecture 
Non-linear analys There are many ases of importance t ndust! where the 
hysical and mathematical formulations for a situation are fairly we established, but 
the purely mathematical theorie nd techniques needea do not exist and consequently 
should be developed The field of non-linear analysis furnishe many current ex 
amples stemming from problems of applied mathematic n the field f dynami 
hydrodynamics, aerodynamics, et I refer to problen n which the usual principle 
of superposition do not hold; problems in which the equations encountered are not 
linear and which when approximated by linear means lead ti olution which are 
incompatible with the underlying physical problems. The development of mathematical 
methods and techniques to cope with such problems is the ibject much research 
these days Improved methods for producing reliable prediction Wi be developed 
in time In the meanwhile the industrial and applied mathematician mu io the best 
he can t ope with such problems and to construct the simplified model I ich a 
way that he can make the most reliable predictions possible 
As a familiar example of a type of problem in this field, I mention a cla of problems 
relating to the motion of ships. In such problems one form ofa differential equation 
which must be studied in order to understand the behavior of shi 

(1) x + (ax + bx?) + cx 0 

where x is the angle of roll of a ship in still water, Whenthis equation is approxi 








mated by an equation of the form 


and solutions of (2) are obtained for certain boundary conditions, it is found that the 
solutions so predicted are completely at variance with the known behavior of a ship 
of the type and conditions under study The replacement of equation (1) by equation 
(2) for purposes of calculation and study oftentimes disrupts the soundness of the 
anchorage so that the solutions deduced are practically worthless in the field of study. 


As another example, I cite the familiar Van der Pol equation 


9 
(3) z+ & ts lypx+xe=0, 
which is encountered in the study of the simple vacuum tube. Theory shows that on 


the slightest provocation the current x must go into a well-defined oscillation. Explicit 


solutions of this equation (3) are often too complicated to be useful or else are not 


known For these reasons the approximation 

(4) x kx + x 0, 
is sometimes used However a study of the solutions of equation (4) indicates that 
the current x will increase indefinitely once it has been started. This, of course, is 
contrary to the facts. Again the trouble lies in an inadequate anchorage of the mathe- 


matical model, equation (4), in the actual field under study 


Use of simulators In the development of such devices as automatic pilots and con- 
trol systems; in the study of the electrical and mechanical stability of interconnected 
systems; and in the study of aerodynamic stability, the typical dynamical situations 
involved often lead to non-linear differential equations of difficult types. Attempts 
to replace these equations by linear systems commonly lead to grossly incorrect 
predictions as tothe behavior of the physical systems under consideration. In order 
to avoid the pitfalls of such poor anchorage, engineers, scientists and mathematicians 
have resorted to the construction of various types of simulators to carry out the de- 
sired solutions. Such methods have been fruitful, but they have also frequently failed 
because of the failure of one or more of the components of the simulator to faithfully 
simulate the desired equations. This is again a failure due to poor anchorage of the 
system in the original setting of the problem Experts in this field all agree that 
when recourse is made to the use of simulators to solve such problems, special care 
must be made to insure suitable anchorage; and such experts go farther by recom- 
mending that an independent dynamic analysis be made of the problems which the 
simulators are designed to solve in order to furnish as far as possible a sound foun- 


dation for and check of the mathematical and physical systems involved 


Bridge design The engineering mathematicians concerned with the design of structures 


ind bridges have over the generations founc much to interest them. To design a safe, 
economical, durable, and useful bridge is an art few men master The history of 
bridge design and construction reveals many tragic instances where the lessons learned 
in one generation have been forgotten in another In every major case of failure 
in modern times, it appears that the mathematical foundation upon which the design 
was generated was faultily anchored in the physics of the problem The Takoma 
Bridge disaster is a recent example Here it appears that the underlying equations 


did not take adequate account of the stability of the bridge under aerodynamic forces 


in relation to its flexibility and other characteristics 





In applying the various theorems and formulas of a branch of mathematics or statistics 


i. 


to an applied field, it is very essential to make sure that all of the mathematical 


premises underlying the formulas of theorems are properly satisfied in the field of 
application Violations of this principle can easily lead to serious errors in the 
primary field of application. It was the violation of this principle which led the older 
classical dynamics and hydrodynamics into a veil of disrepute. The underlying mathe 

matical formulas used sometimes did not take adequate account of such things as 


viscosity, compressibility and the like; and as a consequence predictions made using 
such formulas and related mathematical theorems were sometimes seriously in error 


Serious errors can be made in using the formulas of statistics by not paying strict 


attention to the theoretical premises underlying their validity and use It is very 
often a temptation to use some statistical formula, say, for confidence limit with 
out being sure the theoretical premises underlying the use of such limits are properly 
met in the underlying field The improper use of such formulas may result in a 
completely false picture of the reliability of the data and of the associated predictions 


for the field of study 


The applied mathematician is commonly confronted with a system whose state §S, is 


l 
expressed by functions of observables which are known with varying degrees of proba 

bility It is sometimes necessary to predict the effect of various influences on the 
system and its state. The solution is likely not to result ina unique state, but rather 
in the form of an expected state In other words the solution may be given in the 


form of a state Sp which may be expected to occur with a probability 


When dealing with such situations the applied mathematician may be inclined to feel 


he is working with less rigorous mathematics But ich a point of view can be 
costly, for he must be very careful that the probability formations he uses are based 
on as reliable a theory and data as can be obtained 


nsidered 


4, Nature of pure mathematics in applications. When pure mathematics i 


aS an axiomatic deductive system, the development of the theory is considered to 
proceed from a set P of definitions and postulates--the selection being largely a 
matter of choice--which are not subject to proof within the theory P j tated in 
terms of certain primitive concepts left undefined in the theory The entire theory 
is completely determined when the premises and primitives have been stipulated and 
the principles of logic to be used have been selected Every theorem of the theory 
can then be derived logically from P Pure mathematical proposition ire of the 
form P implies T where P and T involve parameters x y The mathemati 


is not concerned with the truth of P or T, but rather it is concerned with assertion 
of the type “if P is true, then T is true 


This point of view with regard to pure mathematics leads one to raise the question 
as to why mathematics has proved itself of such value in the applied sciences and in 
industry In any field of application the role of mathemati is to render explicit 


propositions Q and assumptions in the field of study which are implied by the premise 


P underlying the treatment of the problem, the propositions Q being hidden and per 
haps heretofore unnoticed. Thus arises the source of the loose statement: “the power 
of mathematics to discover new facts Actually the results obtained are only those 
implied by the premises 
When the terms of a pure mathemati ire given physical definition und the postu 
lates in their physical interpretations can be considered true, then the theorems of 
the mathematics in the light of their physical interpretation can be taken as nece 
sarily true. Such theorems may tell us that if D,, Dg, are done, then S$), So 
will happen 

ae 








The actual function of mathematics in any industrial or applied field is clarifying, 


analytic, but not strictly within itself predictive. The mathematics furnishes a tech- 


nique The predictive power originates in the initial fundamental laws or premises 
assumed for the field under study. If one makes full and proper use of mathematical 
theory and techniques, the failure of a prediction in a field S must be attributed to 


the unsuitability and inadequacy of the premises and observational data involved--that 
is to poor anchorage of the mathematical structure in the field § 


5. The use of the concept of equivalence. The concept of equivalence lies at the 
heart of nearly every physical and mathematical system In the study of physical 
systems, perfect or approximate isomorphic systems play a fundamental role. Equiva- 


lent physical and mathematical models are integral parts of the structures used 


One abstract formulation of the concept of equivalence is as follows: 


Let A, B, C, be elements in the fixed system S being considered (such as abstract 
entities, quantities, ), and suppose a possible relationship that may occur between 
these elements A, B, C, is defined in some manner or other. Suppose that when 


A and B are so related, the relationship is indicated by the symbol — , and the ex- 
pression A£B is written; and that when A and B are not so related, the fact is in- 
dicated by the symbol A&B. Thenthe relationship expressed by the symbol A£B is 
known as the equivalence relation provided that it is defined to satisfy the following 
four | roperties _ 


I. Determination For any pair of elements A and B of §S the relation A&B 
either holds or does not hold 


II. Reflexivity For any A one has AEA 
Ill. Symmetry When A&B, then BEA 


IV ransitivity When AEB and BEC, then AEC 


A is said to be equivalent to B if A&B is an equivalence relation 
Every such definition of equivalence, which is ordinarily not unique for the given 
system, constitutes a division of the elements A, B, C into classes 


All of us are familiar with mathematical examples illustrating equivalence. Thus, in 
euclidian geometry two polygons A and B may be equivalent in the sense that they 
have equal areas, but they may or may not be equivalent in the sense that they have 
the same number of sides or angles Again, any two proper conics A and B are 
equivalent in the sense that one may be transformed into the other by a projective 
transformation; while onthe other hand, A may or may not be equivalent to B inthe 
sense that A may be transformed into B by means of a rotation and a translation 


In general, in geometry or physics, whether two quantities are equivalent may be 
determined by calculating certain numbers called invariants which are associated with 
the object under study 

4 related mathematical concept of considerable importance in its own right, as well 
as in its application, is that known as isomorphic systems 


Consider two systems A and B each consisting of a set of elements and a set of 
operations on these elements. Suppose that each system is closed with respect to a 


given system of operations in the system. The two systems A and B are said to be 





isomorphic, or abstractly identical, with respect to these operations if there exists 
a bi-unique one-to-one correspondence between the elements of A and B, such that 
any formal combination of, or operation on, the elements in A corresponds to the 
analogous construction with the corresponding elements in B 


In any specific study of physical phenomena, as in the applications of mathematics, 
the elements a, b, c, of a system §S are identified with specific physical objects, 
numbers, quantities, entities, or the like, and the definitions of equivalence used must 
be carefully given in terms of these elements and the physical systems to which they 
belong Each definition of equivalence used constitutes a separation of the set of 
elements a, b, c, and associated physical phenomena or systems into classes 
This separation into classes may be of considerable physical significance, each type 
of equivalence often having associated with it an extensive physical theory 


Both in mathematical theories and in physical theories many non-isomorphic types 
of equality have been and can be defined. Each type of equality may have associated 
with it an extensive theory. Thus, in the theory of equivalent linear electrical circuit 

two 2-pole networks may be equivalent in the sense that for all frequencies they have 
identical driving-point admittances (or more generally, for 2N-pole networks, that 
they have identical characteristic coefficient admittance matrices), yet they may not 
be equivalent in the sense that they have the same number of independent mesh circuits; 
or, are structurally the same; or, are equally economical to operate; or, are both 
readily physically realizable; and that if A and B, respectively, are their network 


matrices, A and B may or may not be equivalent in the sense of matric congruence, 


In the design of structures, ships, dams, and many other things, it is necessary to 
make accurate predictions as tothe characteristics, cost and performance of various 
proposed designs Often it is necessary to use models in this type of work The 
results of the studies with the models are then used to predict the performance and 
characteristics of the proposed prototype Great care must be taken in using such 
methods, and the mathematical structures involved must be studied carefully with 


respect to their anchorage, both in the model and in the prototype 


Theoretically, in designing a model test and in constructing the models for the test, 
every physical parameter of impcrtance in the prototype must be considered and taken 
into account It is not sufficient merely to make the dimensions of the model and 
prototype proportional If a theory of similitude is not available, then one must be 
developed. Such a theory must serve as the basis for designing and constructing the 
model from its prototype dimensions. It must serve also as the guiding principle in 
interpreting the measurements made on the model in terms of the corresponding 
measurements as predicted for the prototype. If this theory of similitude is a perfect 
one, and the application of this theory is also perfect, then there is a one-to-one 
bi-unique correspondence between the various physical parameters of the prototype 
and the corresponding ones of the model; and every operation on or with the prototype 
has a corresponding operation on the model In other words, any characteristic of 
the model as predicted through the isomorphism thus stated must be a true character 
istic of the prototype, and vice versa. (This is an example of the equiva.ence property, 
if A&B, then BEA. ) 


In practice it is impossible to apply perfectly ich a theory and expect the results 
predicted through the model test to be completely true for the prototype. A measure 
of this perfection lies in a comparison of the actual characteristi of the prototype 
with those predicted by the model studie Consequently the aim in all such work is 


the construction of as nearly perfect an equivalent system of models as is humanly 
possible In other words, these models must be as nearly isomorphic with their 


corresponding prototypes as possible 








In one type of theory of similitude which is widely used in aero- and hydrodynamics 


a certain set of dimensionless numbers 


are used which must be kept unchanged if true similarity is to exist between the 
characteristics of the prototype and its model. In practice the dimensionless numbers 
nm; of concern are generated from quantities measuring important physical quantities 


; 


‘ 

involved in the study For example, where inertia and gravity forces play an im- 
portant role, the Froude number has been commonly used as one of the dimensionless 
number If true similarity is to exist, every dimensionless parameter 7; referring 


to conditions in the model must have the same numerical value as the corresponding 
parameter for the prototype This commonly means that the model and prototype 
must be completely similar geometrically. From a practical standpoint it is usually 
impossible to preserve the invariance of all of the dimensionless numbers Tj for the 
model and prototype For this reason the experimenters, as well as the designers 
and mathematicians concerning themselves with such problems, must pay unusual 
attention to the nature of the similarities or dissimilarities, and must watch carefully 
the foundation of their mathematical structure in the physics of both the model and 


the prototype 
r wt 


Thus, the use of equivalent systems appears as a fundamental process in almost all 
studies of physical phenomena and in the applications of mathematics to various 
phenomena. Because of the hazards involved in the use of such systems, the industrial 


mathematician must necessarily pay great attention to the tie-in of his mathematical 


tructure with the underlying phenomena under study 


6 Evaluation and analysis of systems In recent years we have witnessed the 
emergence of a new field of applied mathematics, various branches of which are known 
as operations research, evaluation and analysis of systems, and the like. This field 
is capable of very wide industrial application All branches of the armed services 
have made wide use of such field aS a very necessary and useful adjunct to the 
management, planning, and development of large-scale operations and projects and 


as an aid in making important decisions involving large expenditures of money, material 
and human effort. It is in these fields that I have devoted much of my time in recent 
year The foundations of these fields are well established in some respects but not 


in others 


It is only natural that such methods are being increasingly used by management since 
the acts and decisions of executives in a closely-knit industrial organization such as 
we have in this country tend to have very deep and acute influences upon the course 


of events, the flow of production, as well as upon industry as a whole and the indi- 
vidual The resources of this iztry have reached such a state that few industries 
can any longer afford t ise intuition, hunches, and the like, when scientific and 
analysis methods are available which will reduce the risks involved. The margin of 


profit within which industry now must operate is so slim that it is imperative that 


careful analyses be made relating to all important steps 


Most large scale industry has much need for evaluation groups which have as their 


‘ 
inction 


(1) Clarification of objectives 


(2) Clarification of an operational or production process or status 


(3) Analysis of plans, projects, operations, existing methods, and development 





programs with comparative studies and predictions which indicate the di- 
rections in which improvements might be made with existing facilities and 
modes of operation. 


(4) Analyses which indicate the need for exploring new fields and new methods. 


(5) Analyses and predictions which indicate the directions in which investigations 
or projects should be pointed in order to realize the optimum desired results. 


(6) Analyses and predictions which exhibit the courses of action open, together 
with studies of what such actions entail 


(7) Aid in the formulation of policies 
(8) Aid in the technical direction of various tasks, projects and the like. 


In this kind of work which bears upon management decisions, it is extremely im- 
portant that the analyses and evaluations be well anchored in the particular fields to 
which they apply. Where poor anchorage is involved, the price may be extremely 
great tothe industry, which usually cannot afford costly mistakes. Consequently, ex- 
treme care must be made in formulating the foundations for such studies. This 
requires a deep knowledge of the fields involved and the best logical approaches that 
can be devised. The use of incomplete or questionable measures of merit, the over- 
simplification or the over-generalization of a study are apt to be dangerous in this 
field of endeavor 


In spite of the difficulties, no large-scale industry can afford to do without adequate 
scientific evaluations and analyses In this field applied mathematicians, as well as 
engineers and scientists as a whole, can contribute much, particularly when they con- 
centrate diligently on the anchorage of their work in the field of study 

7. Conclusions. To summarize, I have discussed in a broad way some of the prob- 
lems which the applied mathematician must face in industry, and have placed particular 
emphasis onthe great importance which he must attach to the anchorage of his work 
in the underlying fields I have outlined a certain fundamental process common to 
most applications of mathematical methods In this process the use of equivalent 
systems appears as a fundamental process Unless the foundations of such systems 
are firmly rooted in the primary field of study, the results are most apt to fail-- 
without proper anchorage the work of the applied mathematician will fail and perhaps 
lead to costly errors 
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The commands are made to refer to the words 


way. Addresses are assigned to the areas (bins) of 
Let us assume that the 


then be assigned the addresses 0 through 999 


word to indicate the operation (addition, 


remaining twelve digits to specify four of these 


We can now 


It instructs the control to add the number at 
out destroying 
what is already there) 
Thus the threes 


then to go to address w for 


first arguments of the ommand 


operation itself and the last argument with executing a sequence of operations of 


it is a part A good notation for expressing the 


command is 


“[x]” means the contents 
the contents of bins x and y is to replace the prior 
held the 42 97, and 


was executed, 


where 





if bins 1,2, and 3 
oe * Fie 


pectively after the command was 


numper 


mand “+ they would hold 


executed 


By means of the above notation we can easily 


ommands: 


As in the addition command, the last argument 
ft 


the location of the command to be executed next 


will start a problem 


t Thereafter each command 


address 0 


the memory which store 


division, etc.) then there is 


describe the commands. ’* There is one 


address x to the 


either of these numbers) and place the answer at 
have to 


arithmetical portion of the 


of bin x, and the 
contents of bin z.* 


323 respectively 


(when so instructed by the operator) by exe 
it exec 


to be operated on in the following 


words. 


memory can hold 1,000 fourteen-digit words; the bins will 
If we allot two digits of a command 


room in the 


three-digit addresses 


¢ 


for addition which is written 


number in bin y (with- 


address z (erasing 
next command to be executed. 


do with the arithmetical 


£ 


which 
addition 


arrow Signifies that the sum of 


For example, 


before the com 


the numbers 42, 97, and 139 re 


remaining arithmetical 


j a 
Pe} - [y]—t 
[x] Ly} ~(z] 
[x] + fy)—-[] 
r 31/2 

[xt] <L ge [SI 
[x] —[z] 


each command (“w”) specifies 


We will assume that the control 
iting the command 
utes will tell it where to find the 





' We give hers internal 


only the 


into and out of the machine The reason for tl 


we wish to confine our attention to what goes on 
ignore ich factors as decimal points and rangs 
a machine, of course, all of these things are ve 
plexity of the programming The general princip 
out more clearly in the less complex situation we 


9 


~ This is a f 


ommands, omitting those for 

his that, in the 
inside the 
- 


modification of a notation proposed by 


ordering information 


interest of brevity, 
shall also 


actual use of 


machine We 
numbers In the 


; 


important and add to the com- 


f programming, however, stand 
‘al with here 
Hartree, Op. cit., p. 113 














next command The computation may be stopped by using the command 
I 
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orde! the t t modity the ommand at y by ibstituting in its i’th argument 
pla * or 4) the number at addre x, and then to go to z for th next 
omman These are riled g i ommand because they involve the logic i! per 
ation oO nditiona witching ar ibstituting needed to control and organize the 
omputatior 
2 The Log f Programming 
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alle gram sequenc‘: r rogram The activity f onstructing ict equence 
" 1 rogramming. It the maj purpose of thi aper t« eve le prir le 
f programming, that how how programs may be onstructed it of ymmands, 
numbe ina otf progran 
Our task an { n ! 1 by ikKing riou y th fact that a macnine languag‘ " 
inguag It f é ur tificia anguage, xpre y designed rs ur 
ticular ymiputer r mtrast to a anguagt ike f nglish wnicn 1 i natura inguage 
that ha leveloped gradually without much explicit ynstructior Ther ire, never 
the le mportant IY t ul tructure between the machine language and the 
English inguag* Thu ! ler ow we make ommand in English e.g., “Add 
th first twenty number on the boar " Such a ommand is ompounded out of 
ma r nguist part the word ‘ow in a broad sense a program in the ma- 
hine inguag* a ommand or imperative to the machine to perform some mort 
Oo! é ymplicated operation and it also is compounded out of smaller linguisti 
art (numer a and ommand words) Thus instructing a machine is like com 
manding a human being in that we e a language with a limited vocabulary of words 
and expr lifferent commands by ombinations of these 
For th reason our examples have been hosen for their didactic value rather 
than for their realism 





Now a problem to be programmed is generally first formulated in ordinary English 


and mathematical symbolism In the next section, for example, we are given the 

task of programming the imperative “Add all the numbers from a, to ay.” It follows 
s 

that programming is a kind of translating: the function of the programmer is to 


translate a set of instructions or commands from ordinary language into the language 


the machine understands And from this it in turn follows that we will be aided in 
our analysis of programming by any available relevant knowledge concerning the 
tructure of languages Modern ymboli logic provides us with some information 


of this kind, which we shall utilize 





That mathematical logic is relevant to programming will no doubt appear surprising 
to many readers, as perhaps does the use of “logic” in the title of this paper. Ad- 
mittedly, the subject matter of traditional logic (rules for deducing true conclusions 
from true premises) is not relevant here, ince commands are neither true nor false 
and are not deduced from one another as are declarative entence However, con- 
temporary logicians investigating the nature of valid inference have found it necessary 
to study the grammar or syntax of language i.e., how compound symbols are con- 
tructed out of their linguist parts (primitive symbols) As we have just seen, a 
rogram is a combination of word (in the machine language) which constitutes a 
ompound imperative Hence our problem is in important respects similar to one of 
the problems the mathematical logician studi It is therefore natural that some of 
the oncept employed in logi (e.g., propositional function, free variable, bound 
variable, and quantifier) are useful in the theory of programming 
We must emphasize that there an important difference between machine languages 
und the artificial language isually investigated in symboli logis The latter are 
arative languages,+ whereas machine languags contain imperatives (commands); 
it is with commands that programming is mainly concerned However, as 
we have a‘ready indicated, we are not presenting a logic of imperatives in the ensé 
of providing ruls of deduction for mmand Rather taking advantage of the 
imilarity of the grammar of imperative to the grammar of declaratives -- we ar‘ 
applying gica oncept jeveloped f« é irative in our analysis of imperatives 
2 The Structure f a Simple Program 
In thi ection we will analyze a imple program, and in the subsequent sections 
how how ich imple programs may be ompounded to give mor‘ omplex one 
Suppose we wish the machine to perform the simple summation 
I 
(1) , aa 
What we want the machine to do is expressed in ordinary language by the command 
“Add together the number a; through ay.” Since, however, there is no (single) 
ommand in the machine language for repeated summation, we must make the im- 
perative more explicit before we an translate it into this languages 
One way of forming the required sum is to begin with the quantity 0 and successively 
add the a’s to it Such a proce i alled a recursion or induction and is best 
Variou logi of imperativs have been proposed, but the writer did not find 
them o any i1¢ n the problem at ha 
2Q 
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functions For example, “+ (18,6,3,139)” can be expressed in the machine language 


as a fourteen-digit number, but “+ (x, y,z,w)” cannot since the machine language 
described in Section 1 contains no variable symbols.‘ We did employ variables to 


explain the machine language, since we wished to describe the action of the computer 


for any definite command of a given form For the same reason command functions 
and ordinary numerical functions as well (e.g., “g,;" in bin 1 below) are essential 
in programming Hence we will extend the concept of machine language to include 
such expressions This will cause no confusion as long as it is remembered that 
before a program is actually given to the machine (i.e., before the task of trans- 
lation is fully completed) the variable symbols should be removed (see Program II 


To program (c), then, we must write a command function in the machine language 


which is a translation of the English command function “comput gj by adding gj_) 


i 
to aj.” Suppose that the a;'s are stored in the bins L + 1, L + 2, iae we E Oe 
that a; is stored at L + i) and g; is stored in bin 1, Then the following is a trans- 


lation of (c) into the machine language: 





g 
2 + (1 L. + i ] ) J y + a 
: _— 5 j Si-1 1 
We may call this part of the program the computation part The symbols on the 
left describe (by means of the variables) the contents of the bins at each stage of 


the computation Note that the variable “i” i hosen so as to describe the contents 
of the bins at the end of the execution of this part of the program rather than at 
the beginning The notation on the right is to help the programmer keep account of 
the role the various words play in the computation 


Because the control understands and executes on!y definite commands, not command 


functions, we must instruct it to specify properly all command functions into definite 
ommands before it reaches them This, in fact, is what is commanded by (a). 
(a) tells the control to keep account of i, the number of times the command function 
at address 2 is executed, and for each value of i to make the proper substitution 
into this command function A translation of (a) into the machine language is thus: 3 

3 l 

+ A Constant 

2 L IL. + i is the location of a, 

f Temporary storag« 





8 +(3, 5,6, 9) f L. + i 

) S(6,2, 2, 2) 
* A language containing variable symbols and which the machine can in a sense 
understand is suggested in the last paragraph of Section 5 
- : 

The situation may be described more precisely in this way: the programmer uses 

a metalanguage containing variables to talk about an object language ( the language 
understood by the machine) which do not contain variable It is unnecessary for 


us to distinguish these two languages in our subsequent exposition, however. 


























Since the functior tk part the rogram is to specify properly the ommand 
func tis I mputation part int finite commands, we will call it the specifi- 
ation part tk gram rt " ymmand is writter o that the ontrol will 
witch t tr 8 itatior art I UF rogram wher t na ompieted tnis art 
vt nav ot tf rT arit i) t a’) Ir iOg1 “for very x” A a 
a qua t r t tat tr jJuantit I i tn instan ail [Or wnicn 
th ; tional f t wing it x a man then x is mortal”) holds true 
Simila " ict th mact t xecut the ibsequent mmand function 
f te imbe: tin l the purpose of (b) to stop th omputation 
t tf I . [ f t r on | tional ne tive “if tr n . 
toget! with tf ter I I t ng the computati t A trar ation of (b) int« 
t mact inguag hu r tt nditiona hift ont mmand and I: 
( ¥ f < Itt ntr to 7 
| | ther ty ; a 
Y t j w ( t know yet where tne ontrol nouida g¢ irter 
tn n tal ru r ] Ve w i if part of tr rogram tne¢ I ional 
art t inctior to nti tk repetitior f the e in accordar with 
the g ‘ r ‘ 7 
I bta i ng rogram which 1 1 translation of (2) into the machine language 
we must omplete the terconnect of the three parts by s} f last ar- 
gument of the mmand at 2 Su ne onditional par yr =the I hould be 
executed afte tk mputation part ommand 2 should read 
2 + L, + . 
The sequence of word n the bi through and L + through L + I nstitutes 
the program for the problem specified We hall number it Program | 
Actually Program I do no onstitut 1 definite set of instructior to the machin 
Thi the ase becaus«¢ ich quantiti¢ i gg, L, and th i” have been left un 
pec ified In the termir gy used irlier Program I is a command function, not 
a definite ommand (Cf. the or imperative (2) in thi respect.) Since it ( 
em best to r rve the word ommand” for the imp] ommands introduced in 
Sectior we wi ay that Program I i 1 program function and any complete specifi 
ition of it a definite program “Program” i used indifferently to refer to 
tner I nstr ct a ma hine to olv 1 problem we must, o! ourse, give it a 
lefinite program ince only definite programs (not program functions) are under- 
tood bv th ntr ; 
We off as an example f a definite program the f S ification « Program I, 
which w ha i Program II 
t Initial value of g 
2. +(1 , 11) g = g; + a; 
3 Phi equence of words dot no juite onstitute a complete translation of (a) 
net t dos not fix the initial valu of i and g This is done by specifying the 
ontent of bir l and 3 to be O before the omputation commanded by (2) starts 





3 0 Initial value of i 


4 ] Constant 
5 1] I 
6 Temporary storage, initial contents 


irrelevant 


“ 
+ 

oe 

Ww 


4,3, 8) i —»ji-] 


12 123 ay 
13 247 ao 
14. 181 a. 

3 
15. 49 a, 


If we add to this sequence of words commands for getting to it from it-- 


0 ey My Pa 
~~. F 


--and express all of the words as fourteen-digit numbers in the machine language, 
we will have a complete problem ready to give to the machine, The machine will 
start with the command at 0 and then switch to the command at 7.! It will cycle 
the commands 7, 8, 9, 2, 11 in that order four times and then go to 16 and stop, 
signalling the operator that the problem is finished At this time address 1 will 
hold 600, the answer to the problem The reader may find it helpful actually to 
follow the control around the feedback loop a few times, noting how the contents of 
bins 1, 2, and 3 vary during the computation 


It is important to notice what happened to command 2 in the specification of Pro- 
gram I into a definite program A blank is written in the second argument place 
since no relevant number could be placed there: no matter what number is placed 
there the program will substitute the correct value of *L + i” in its stead at the 
proper time. Thus in converting Program I into a definite program we did not assign, 
and could not significantly assign, values to these occurrences of the variables *L” 
and “i”. Since the specification part of the program does this, Program II is a defi- 
nite program despite our failure to specify these variables (Cf. in this regard the 
fact that “For every x, if x is a man then x is mortal” is a definite proposition 
even though it contains three occurrences of the variable “x”.) In contrast, the 
occurrence of the variable “L” at 5 must be specified to convert Program I into a 





1. Note that the conditional switch of control command becomes an_e unconditional 
switch of control whenever all the arguments are identified There are other ways 
of producing the same result, e.g., in this case “C(3,3, 8,7)’ and “C(5,4,7,7).” 
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ie terminology of symbolic logic we can say that the 








"7 and “i at 2 are bound occurrenc: in Program I 
iriable [ and “I” at 5 and 10 are free o ir 

ordinary mathematical distinction betweer variables 

e of a variable is bound or free is relative t i 
Thus though the occurrence of i” at 2 is bound rela 


relative to the omputational part of this program; 


. ) 
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the control through the program from the entry point 7 to the exit point @. A double 
arrow is used in the specification box to symbolize the fact that it is the function 
of this part of the program not only to advance i by 1 for each cycle but also to 
specify all occurrences of free variables in the computation part of the program that 
must be respecified for each new value of i (in this case “i” and “L” in 2), (See 
2 on page 44. ) 


footnote 


4 Integrating a Trajectory 





Let us now program a more complicated problem to show how the principles estab- 
lished in Section 3 may be used The problem the first general purpose electronic 
computer (the ENIAC) was built primarily to solve was that of the trajectory, and 
it will serve as a good example 





We begin with the differential equations for the path of a projectile in motion:! 
=f -Ex! 
yi -Ey' -g 
where 
“RY G = 2 
i. 3 G(s) : s “"y (x') + (y')*, 
C 


x and y are the horizontal and vertical positions, a prime denotes differentiation 
with respect to the independent variable (t), g and h are fixed constants, C is a 
constant for a given shell, and G(s), the ballistic drag function, expresses the re- 
Sistance of the air to the shell as a function of its speed. These may be reduced 
to the following set of simultaneous first order equations: 





u x* 

y y! 
u! -Eu 
v' = -Ev - g 

E ety G( Vu2 + vy") 


To get an approximate solution to these equations by digital computation we replace 


them by corresponding difference equations. These are, symbolizing e “hy by “F(hy) 
and the time interval of integration by “A”: 





2 2 1/2 
~n-1 (uy, _} tv, 24) 
E F(hy,, » G(s.- ) 
C 





For the derivation of these equations see G. A. Bliss, Mathematics for Exterior 
Ballistics, John Wiley and Sons, New York, 1944, Ch. 2 
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‘ ‘ +A 
xX Xy + Un.;4 
y, y, + Vyn_14 
J ] I L..7@ 
V Vv, I v. ,4-g 
It wi be noted that these equations are recursively formulated, so that given the 
alues of the appropriate quantities at the beginning of an interval (i.e., tp-1, Xn-1, 
Yr — and vn-1) we can by a series of arithmetical operations produce the 
orresponding quantities for the end of the interval (i.e., th» Xn» Yr» Un> and Vn) 
Ther ire two function involved in the calculation The drag function G(s) is de- 
termined empirically and hence is provided in the form of a function table The 
exponential function F(hy) is normally produced by solving a differential equation as 
the omputation roceeds, but in order to illustrate how the same program function 
in be u 1 more than on in a single cycle of computation we shall assume that 
t too is given in the form of a function table We then need a program to look up 
value of y given function f(z) in a table, so let us write that next 
We wi i ime that values of f(z) are given for successive integral values of z be- 
ginning with O and are stored in sequence beginning at address M. If we assume 
further tha in integer’ then f(z) is located at M + z The following program 
inction (Program III) will look up th value of f(z) in the given table: 
2 Argument 
22. M Location of f(0) 
23. f(z) Function value to be found 
24 +(2 22,21,25) 21 M + z 
c ’ 
25. S(21, 2¢ 2t 
vA r'(M + , 223 6B) Transfer f(z) to 23 and nitt ontrol to B. 

It rep! nted by the low diagram of Figure 2, which has a specification box 
and computation box in that orde 

We in now ogram the trajector\. A tr equations show, the first tep i to de- 
termine F(h ) and G(s, He we must ymmand the machine first to specify 
Program III into a definite rogram to look up F(hy,_,) The control should then 

’ 
We make t A imptions to simplify the program even though they are not real- 

ti For th ime reason we d 1ot interpolate and we use only first order differ- 

n quations (3) instead of on of a higher order In actual practi efficiency 
with gard to mputation time and memory space would require us to do all these 
thing But iding them would add to the omplexity of the program without in- 
troducing any new principl« of coding, and as we pointed out before it is more im- 
portant for o1 xample to illustrate basi 


principles 


clearly than to be r 


realisti 





switch to Program III and after finding F(hy,,_}) it should switch back to the main 
routine to specify Program III into a definite program a second time, this time to 
look up G(sn-}1) After passing through Program III again the machine is finally 
ready to compute ty, Xn, Yn. Un, and vy This whole process should be repeated 
until some criterion is satisfied. Suppose we wish to continue the integration until 
the shell “hits the ground,” i.e., until y becomes negative (or zero). Then a con- 
ditional command to control the cycling will conclude the program. The structure 
of the whole program (Program IV) is represented by the flow diagram of Figure 3. 
Note that, as this diagram makes clear, Program III is included in Program IV as 
apart. £8, and Bo, the two points to which the control returns after executing Pro- 
gram III, will be determined as the program is written. 


By using the flow diagram we can write the program in a straightforward manner. 
Let us assume that 100 values of F(z) are stored at addresses 100 through 199. 
The occurrences of free variables in Program III to be specified to make this into 
a definite program are “z” at 21, “M” at 22, and “gf” at 26. (“M” and “z” in 26 are 
bound and hence need not be specified. ) 


31. h 

32. yn This bin will hold y,_, at the beginning 
of this step of integration, and y, at the end. 

33. 100 Value of M for F(z) 

34. By Location of command to which control re- 
turns after looking up F(hy,,_;) 

35. X31, 38, 21. 36) oo ee in 21 

36. T(33, 22,37) 100——> M in 22 

37. S(34, 26,3, 24) B, ——~8 in 26 


The last command switches the control to Program III. Since the first command 
is located at 35, x 35 


This same process is to be repeated for G(s,_,), after F(hy,_,) has been transferred 
to temporary storage. Let us assume that 100 values of G(z) are stored in bins 200 
through 299 


38. up This will be u,_; at the beginning of the 
step of integration. 

39. *. This will be v,_,; at the beginning of the 
step of integration. 

40. Temporary storage 

41. Temporary storage 

42. 200 Value of M for G(z) 

43. Bo Location of command to which controi re- 
turns after finding G(s,_1) 

2 
44. X(38, 38, 40, 45) 40: u,-] 
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9 
45. X(39, 39,41, 46) 41: Yn-] 

2 
4° +(40, 41, 41, 47) 41 s,- 
47. Vv (41, 21, 48) $4.) es, in 21 
48 1T(23, 40, 49) 40: F(hy_,) 
49. T(42, 22, 50) 200 ——.M in 22 


50. S(43, 26, 3, 24) Bo —+ 8 in 26 


Note that 8B, 44 
The execution of the second specification of Program III will leave G(s,,_;) in 23. 
We can now proceed to compute t, x, y, un? and Vv,» using the difference e- 
n i r 
quations already given: 
51 th This will be ty-1 at the beginning of the 
Cc y¢ f 
59 This i7] . } > . , f , 
02. X, This will be x,_, at the beginning of the 
cycle 
53. A 
94. ( 
95. g 
Df ] & 
57 Vv, ya " 
58 +(5 53,51, 59) t t ; 
ni-—— n-] 
59. X(38, 53, 56, 60) 5€ u~—« 
60. X(39, 53,57, 61) 97 t. , 4 
6 +(52, 56, 52, 62) a 
-9 s > 6&7 ) - 
62 (32,57, 3 63) y y, 
63. X(40, 23, 40, 64) 40 F (hy, ;) G(s,_)) 
Yn-1 ~ 
64 (40, 54, 40, 65) 40 E. 
65. X(40, 56, 41, 66) 4) E u ray 
6 E (38, 41, 38, 67) } u 
_— ht n- 
67. X(40, 57,41, 68) 41: EE. hen ra\ 
68. X(55, 53, 40,69) 40: gd 
69. +(40, 41 1, 70) 41 E j A+gS 
40,41,4 = os g 
70 (39, 41,39, 72) = j 
= a. B- 
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This completes the part of the programming represented by the computation box of 
Figure 3. Note that Bo - 58 


Finally, we code the conditional box: 


71. O Constant 
72. C(71, 32,35, 3 ) If 0 < y, then control to 35 
If 0 2 y, then control to & 


This part of the program illustrates an interesting point. The path of the control 
is determined here by one of the answers to the problem (y,) rather than by a num- 
ber introduced primarily for purposes of control (e.g., the control number i in Pro- 
gram I). Thus there is a feedback from the numbers being processed to the com- 
mands. A like feedback occurs in Program III: what definite command the command 
function at 26 becomes depends on the data In general, the path the control follows 
through the memory depends on the numbers being processed as well as the com- 
mands and control numbers which instruct and guide the computation It is for this 


“ 


reason that we defined “program” to include words expressing data and answers as 


well as words (commands and control numbers) whose function is to direct the com- 
putation A program is thus a sequence of words which contains instructions for 


its own transformation 


Program IV does not constitute a complete formulation of the exterior ballistics 
problem. Apart from the fact that we oversimplified the process of determining the 
functions G(s) and F(hy) and used only a first order integration technique, Program 
IV is only a program function, not a definite program. In actual practise it would 
be employed repeatedly for different initial conditions (muzzle velocities, angles of 


fire, etc.) and different shells Hence it would need to be embedded in a program 
of wider scope which would bind the free occurrences of its variables: “up,” “vo,” 
“2° © 3.* on This program would also need to provide for reading in data and 
writing out the answers, as well as for transforming the trajectories computed into 
firing tables, the end result of the computation Though it is beyond the limits of 


this paper to code the complete program, we shall make a few remarks about the 
construction of such programs 


5 Program Composition 





Our discussion of how best to construct complicated programs will be facilitated if 
we can refer to a further example Suppose we wish to instruct the machine to 
compute 
J / I \ 
| | ae 
(4) \ ca; I 


We can most simply code this problem by using Program I to compute the quantity 
in parantheses and adding to it a program corresponding to the multiplication sym- 
bol “tT” As before, the problem is best formulated in inductive terms Define the 
function h; recursively by 








This example is interesting because it involves two quantifiers and hence its com- 


putation involves a double recursion Since, in the formula to be computed, the 
multiplication symbol “tT” is of wider scope than the summation symbol “2” we 
need a program with a feedback loop which will enclose the feedback loop of Pro- 
graml. The flow diagram of the required program is shown in Figure 4. Program 
I can be made a part of this program by specifying @ in the command at 11 to be 
the location of the first command (in order of execution) of the computation box of 
Figure 4. The specification box of Figure 4 is the only one that needs special com- 
ment 


In general, the function of a specification box is to bind the free occurrences of cer- 


tain variables in a particular sequence of words (Similarly, the function of a quanti- 
fier in logic is to bind the free occurrences of a variable in a particular sequence of 
ymbols.) In this ase these are the free occurrences of the variable “I” and of 
the other variables (“L” and the initial values of “i” and “g;") which need to be re- 
specified for each new value of “I” in the program between the point marked “7” 
and the end of the onditional box in Figure 4 Thus the specification part of the 
program must keep a count of the major cycling (i.e., of I) and for each major 
ycle specify the occurrences of “L” (note that the location number L is now a 
function of I) in 5 and the initial values of “°i” and “gj” in bins 3 and 1 If any of 
these variables were free in the computation and conditional boxes of Figure 4 they 
would have to be pecified there also, as would be the case if, for example, the 


omparison made by the conditional box involved a quantity whose location was a 
inction of I 


This last example and Program IV illustrate a general method of constructing pro- 
grams: that of coding program functions and later combining these to give larger 
program We in take the fundamental element in this construction to be the pattern 
hown in Figure 4 Other patterns may be regarded as degenerate cases of patterns 


like this (e.g., Figure 2 lacks the conditional box and the feedback loop) or as com- 
r oth 


f lem in which the elements are hooked together at the exit and entry 
points (e.g., in Figure 4 the “7” in the exit circle following the specification box 
and the “q@” in the entry circle preceding the computation box connect the diagram 


Let us consider in more detail this method of programming Most problems solved 
on a computer an be broken down into a number of sub-problems. For example, 
the computation of firing tables involves interpolating, solving difference equations, 
using a function table, printing, et Program functions can be written for these 
fundamental processes and then used repeatedly. Since these programs will be lo- 


cated in different parts of the internal memory for different problems, most of the 


addresses must be written as free variables When a specific problem is to be 
olved it is not programmed as if it were unique but is dealt with in the following 
way First, the program functions corresponding to the fundamental processes in- 
volved in solving the problem are selected from a library of programs. Second, 


a statement is made in the machine language of the unique features of the problem; 


e.g., a description of a particular set of equations is given and the initial conditions 








listed. Finally, a combining program’ which will bind all of the free variables (in- 





cluding the addresses) of the sub-programs is added to these programs to give the 
desired definite program 


There is a very practical reason for employing this method of instructing a com- 


puter An electronic digital computer is extremely fast; for example, the ENIAC 
can compute a trajectory faster than the shell itself can fly. It is important that 
the increased speed of computation provided by electronic machines be matched by 
a corresponding increase in the speed of programming. This is possible if pro- 
grams are constructed by combining already prepared program functions. It should 


be noted that the specification of program functions into definite programs is not 
accomplished by rewriting the programs by hand outside of the machine, but is done 
by the machine under the direction of the combining program. This is significant, 
for since most of the addresses in the programs are represented by free variables 
the task of specification is a substantial one Thus a large part of the work of pro- 
gram composition is done by the machine 


We shali conclude this paper by briefly suggesting a way of using the machine to an 
even greater extent in program composition What we have said about programming 
can be summarized as follows The programmer begins with a formulation of his 
problem in the ordinary language of mathematics and English (e.g., (1), (2), (3), 
(4)). He then breaks down the given statement into elementary parts and translates 
each of these parts into the machine language; for example, Program I is a trans- 
lation of a part of expression (4) Finally, he writes a combining program to direct 
the machine to combine the resultant program functions in a way determined by the 
arrangement of the parts of the original statement of the problem; for example, since 
the “tT” of (4) is outside of (1) the programmer makes Figure 4 enclose Figure 


Under this procedure each program synthesis is treated by the programmer as a 
unique problem. Now it seems possible to instruct the machine by means of a single 
translation program to synthesize many definite programs under the direction of the 





original statements of the problems. Suppose we constructed a programmer's language 


in which expressions like (1) and (4) could be written Such a language would be 
much closer in structure and symbolism to the ordinary language of mathematics 
and English than the machine language described in Section 1, and hence easier for 
the programmer to use in stating problems to the machine. The machine could be 
taught to receive instructions in this language in the following way The symbols of 
this language could be represented in the machine languag: The programmer would 


translate simple expressions in this language (such as (1)) into program functions in 
the machine languag: Finally, he would prepare a translation program (in the ma- 
hine language) to instruct the machine to combine these program functions under 
the direction of expressions like (4) to make definite programs Thus we can teach 
the machine to make the translation from the programmer’s language to the machine 
language, and once this is done the more convenient programmer’s language can be 
used for instructing the machine to solve problems 








See Goldstine and von Neumann, Ibid., Vol. III. 
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ANALOGUE COMPUTATION 





Arvid W. Jacobson 


The subject for the February (1950) meeting of the Society was Analogue Computation 
and its Engineering Applications. The meeting was followed by four conference 
sessions at which Professor S. H. Caldwell, Director of the Center of Analysis, 
M.I.T., led the discussions. The meeting also served as an occasion for the opening 
of the Wayne University Computation Laboratory. (Wayne University had earlier re- 
ceived as gifts from M.I.T. the original differential analyzer built by Dr. Vanevar 
Bush and the cinema integraph.) A total of 500 local engineers and research people 
attended this series of meetings 


The conferences began with the description of the differential analyzer, and continued 
with the discussion of its application to the solution of problems in mechanical vi- 
brations, in electrical transients and servomechanisms and concluded with a session 
on scientific problems. In this report of the conferences it is not feasible to repro- 
duce the discussions in full. However, we shall briefly describe the differential ana- 
lyzer, show how problems are set up on it and solve one or two engineering prob- 
lems. (See bibliography references at the end of the paper for more detailed dis- 
cussion, ) 


DESCRIPTION OF THE DIFFERENTIAL ANALYZER 





The differential analyzer is a machine for solving ordinary differential equations. 
In it the value of a given variable is represented by the amount of rotation of a shaft. 
It thus operates on the physical quantity of angular rotation, and hence is called an 
analogue machine. The shafts are so interconnected that they rotate in accordance 
with the functional relations expressed in the differential equation. 


The operation involved in the solution of differential equations is integration. The 
basic component of the machine is a unit which mechanically performs this operation, 
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' Dr. Arvid W. Jacobson is Assistant Professor of Mathematics at Wayne University 
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across the table, a crank is operated manually to keep the index on the curve. Thus 
the value of the function P(x), represented by the shaft labelled P(x), is continuously 
fed into the machine by the operator. The fact that discontinuous and non-linear 
functions or any experimental data which can be plotted can be introduced into the 
machine enables it to solve problems for which the classical method would be im- 
practical. 


HOW PROBLEMS ARE SET UP ON THE MACHINE 





Let us set up the problem dy/dx = f(x), y(0) = a on the machine. The independent vari- 
able shaft x (“x” represents the amount of rotation of this shaft) is driven by an 
electric motor. The shaft f(x) is turned manually at the input table. Hence we have 
two inputs °x” and f(x) into the integrator and the output “y” is then the of f(x)dx. 
The initial condition y(0) = a is introduced as an initial setting of the integrator. 
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(Dots denote cross gear boxes) 
Figure 5 
The solution represented by the rotation of shaft “y” may be taken to an output table 
by means of cross shafts and plotted, and it may also be connected to a counter and 
printed. 


Another example. 
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A motor turns the independent variable shaft “x 


Since “y” is the solution of the 
equation, and hence initially unknown, we anticipate the solution and therefore the 


motion “y” Thus we obtain the motions f(x) and g(y) from the input tables, which 
in turn drive the adder whose output is d2y/dx2. The output of integrator I is dy/dx, 
and that of integrator II, is “y”. Hence the motion which was anticipated is actually 
provided by integrator II This is a mechanical feedback and the second integrator 
“closes the loop” This feedback is also a machine equivalence of equality of the 


two members of the equation and constrains the variable quantities involved to assume 
values in accordance with the functional relations expressed by the equation. 


In a machine layout of a problem it is essential that one and only one source of 
torque is applied to each shaft In the above example the shaft “x” is driven by a 
motor, The torque to shaft h(x) is provided by an operator, similarly for shaft g(y). 
The adder combines these two torques into one and supplies it to the shaft d2y/dx?. 
Next, the output of integrator I provides (through a motor driven torque amplifier as 
explained earlier) the torque to shaft dy/dx. Finally, integrator II supplies the antici- 


pated torque to shaft “y 


In a machine solution of a problem it is necessary to keep the displacements of inte- 
grators and function units within design limits. Likewise good precision requires 
that the number of revolutions of each shaft be fairly large, while the solution time 
should be as low as possible These conditions are satisfied by proper determination 
of scale factors, which we shall illustrate by a simple example. 


The notation Ax on a shaft means that the shaft makes A revolutions per unit of “x”. 
Hence if A is 256 and “x” varies by 10 units, the shaft will make 2560 revolutions. 


Let us calculate the scale factors in generating e* by solving the differential equation 
dy/dx y Let the scale factor of the independent variable be A, and that of the 








’ 

function e*, B The gear ratio between the two ex shafts is “n”. Then 

; ‘ 

(1) Be* == (Bne*) d(Ax), 
where the fraction 1/32 is the integrator constant “k”. If “x” varies by 5 units, and 
if the solution time is set at 10 minutes while the maximum speed of the independent 
variable motor is about 600 rpm, we get an approximation for A: 5A/600 10. We 
choose A 1024 From equation [1) we get Be* ABnex/32, so that An = 32, or ' 
1024 n 32 and n 1/32 
The limits of the integrator displacement are +40 turns of the leadscrew Hence the 
x , , = x ; 
integrand Bne~ must satisfy the following inequality: Bne~ = 40. With the above 
values for “n” and “x”, we get a maximum value for B, i.e., B = 40 x 32/e°, 
The solution on the differential analyzer is generated in a closed loop. Therefore, 
the error in the solution is also generated in a closed loop and has the same form 
as the solution itself The machine is equipped with units to correct for backlash, 
known as frontlash units The differential analyzer, under proper conditions, yields 
results accurate to one part in 1000 
' 
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ELECTRICAL CIRCUIT PROBLEM 





The circuit is shown in Figure 8, where R is a non-linear function of the current “i”. 


The equation for this circuit is 
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If S is closed at t = 0, we are to find “i” 


as a function of time. The machine as- 
sembly is given below. 
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SUSPENSION PROBLEM 




















In order to keep the discussion from becoming too involved we shall consider a vi- 
brating system with only two masses. The arrangement of masses, springs and 
damping units of the physical problem is shown in Figure 9, 


where x,(t) represents 
the exciting force. The equations of motion are 
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Dividing 
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by the mass and integrating once, we get 
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CONCLUSION 





In recent years new methods in the art of solving problems on the differential ana- 
lyzer have been developed. Notable among these is the method of inverse integration 
without the use of servomechanisms developed by J. G. L. Michel. (See bibliography 
reference at the end of the paper.) The application of the differential analyzer to the 
solution of partial differential equations with special boundary conditions has been 
made by D. R. Hartree. 


The differential analyzer is a powerful tool in the solution of differential equations, 
It can, of course, be used in problems for which the solution can be obtained by 
analytical methods. However, if a large number of solutions are desired for different 
values of the parameters involved, the machine method is the only practical one, 
The same is true if the equation if non-linear or if experimental data are to be used. 


The differential analyzer, since its installation at Wayne University, has been used 
to solve a number of important problems dealing with the automobile, Likewise, 
its application to the determination of the protein fractions of the blood to ascertain 
the presence of disease has been noteworthy. 


Its proper utilization for training purposes, for research work in the University and 
by local industry in the solution of basic engineering problems can materially advance 
the methods of research and bring about fundamental improvements in design and 
greater economy and better quality of products. 
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FLUID FLOW MAPPERS* 





A. D. Moore 


The fields in electromagnetic, heat flow, chemical diffusion, and other problems may 
be simulated by the flow of fluids in suitable channels Despite the pioneering work 
of Hele-Shaw in 1901'!) little work has been done in this field up to now, because of 
the experimental difficulties of his method However, if modern dental stone is cast 
in the form required, using a sheet of plate glass as a mold, it is possible to obtain 
the uniform flow passage essential for accurate work (2) 


The major components of a fluid flow mapper are a plaster slab and a piece of plate 
glass, parallel to the slab and spaced upward from it by a small constant distance 
This creates the necessary two dimensional flow space The combination is immersed 
in a tray of water, and small tanks connected by rubber tubing to holes in the slab 
act as sources and sinks to create the flow The flow is made visible by crystals 
of potassium permanganate and barriers are added as needed to represent the boundary 
onditions of the field 

After the lab ha been made by the methods detailed in reference 2, the fields are 
made visible by flowing water over a sprinkling of colored crystals. The flow ‘patterns 
may be ecorded photographically and are then available for detailed computation of 


potential lines from the streamlines recorded. These pictures are intrinsically inter- 


ting as pattert and color photographs of great beauty have been made 


Distributed sources of flows represent the magnetic field in a region carrying current 
Py 


1 
und may be simulated by a “sand bed”(%)(*) This technique can be used also to 
represent the lution for the torsion of a prismatic beam and the flow of heat 


in a region where heat is generated uniformly 


The equipment has been used to solve a variety of practical problems, as a teaching 
uid, and to verify some theoretical propositions relative to electromagnetic fields (4) 
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A LOCATION ESTIMATOR FOR NON-GAUSSIAN DISTRIBUTION , 














D. R. Elliott 
Abstract 
An expression has been found for an estimator of the location parameter of a non- 
Gaussian distribution. The estimator (g) satisfies the conditions (1) g(X ,,X2,...Xp)= 
g(X ,+U, Xa+t Xnt+U) -U and (2) variance of (g) minimum where X,...Xp, are 
the (n) sample variates ordered according to magnitude and not, in general, just any 
arbitrary ordering, and U is a parameter of locatior Both the median and mean 
satisfy equation (1) but not necessarily equation (2) (g) for a Gaussian distribution 
gives g the mean; (g) for a rectangular distribution give g (smallest variate plus 
largest variate)/2; more complicated expressions are obtained for Cauchy distribution. 
For some listributions, the is the above estimator may be considerably more 
efficient than the average of th ample 
Sometimes in product control work, o1 knows from past experience the form of the 
‘ listribution of a statistica ariabl yut wishes to determine by drawing a sample 
of (n) variates, whether the distribution ha hanged its location (shifted to higher or 
lower values) or is the am It istomary in such uses to determine the aver 
age of the (n) sample variates and t lot the average on a chart along with previ- 
ously obtained values, using time as the abscissa If the distribution does not change 
n time, th ] Sive points on the hart fluctuate about a horizontal straight line. 
\ If the distribution shifts steadily with time then th Icce ive points fluctuate about 
i straight line having a slop Si the ba reason for taking a sample is to de- 
termine whether the distribution ha hifted, it desirable that the fluctuations of 
the points on the hart be ma in order to detect any shift as soon as possible. 
If one s limited in the numb iriat it h disposal, he may desire to us«¢ 
a more accurate method of determining the locatior f the distribution than to form 
the average ven if the a llatior more time onsuming The method of calcu- 
ation presented here may not always be the most efficient method but is more ac- 
} urate thar taking the averagt except n those cas that it reduces to taking the 
average 
t seems obvious that any function (g) of the variat« ipable of detecting a shift in 
the distribution should satisfy the onditior 
g(X,,%o,...%,) = g(&)+U, XotU,...X,+U) -t (1) 
Equation (1) states that if th listributior hifts by an amount U, then the average 
value of an estimator (g) must a increase by an amount | One solution of e- 
quation (1) is 
} g(X1,Xo, Xp) g(X,+U, Xort X,,+U) -I (2) 
Although equation (2) is not the m t general solution of equation (1), it 1 Sulficient 
ly general to include the average ar n mar ases, estimators which are conside: 
ably more accurate than the averags 
’ As a gauge of accuracy, we use the ndition that of all possible solutions of equa- 
tion (2), the solution having th ist iriance is the best Hence, if (p) represents 
“A paper pr! nted at the joint $101 t \merican Mathematical Society and the 
Industrial Mathemat Society at East Lansing on February 25, 1950. Mr. Elliott 
s with the United State Rubber ( mpal Detr it, Michigan 








xX 


joint distribution function of the variates XX. ee a 


the 


fPven. dx? dxn a minimum 
where the integration extends over all values of X,,%2,...%p, 


outside a finite range 


Equation (3) states that the variance of (g) is to be a minimum. 


solution of equations (2) and (3) is easily obtained (see Appendix) and is 











The 
x. x | S xp(x+x, X+X9,...X+Xn)dx 4) 
- 
B81, %2,.--Xn) = ~ 7 p(X+X1, X+Xg,.. . XtXp dx *) 
[t not necs ury that the joint distribution function p(X, Xo,.. .X,) be expressable 
in the form 
@ (X ) P(X2) $(Xn) (5) 
However, if (p) i 1 symmetric function of the variables as in equation (5), (g) is 
ilso a symmet! function Moreover, if (p) is given explicitly with a parameter 
f location @ in the form 
, 
’ »(X 6, Xo-6, X79) (6) 
then (g) in unDila 1 estimate of the parameter 6 
If the distribution is Gaussian, then 
2 2 2 
at 2(X1 + X2 + Xn) | 
ind 
a _ 2 2 | 2 
/ 2([x+x,} [x+x2] +... [K+Xn]*) 
x dx 
g a) 
ff & “ >a 4 ' 
2(fx+x ; [x+xJ + i+x,] ) 
dx 
< 2 9 
2(nX"+2X1X1+X2+ Xn} +X1"+Xo"+ Xn J 
Xe ix 
’ uo 
x ee ——$$<$__ 
r= 5 : _— 2 
/ 2(nX“+2X1X) +X 2+ Xn | +x +X» + Xn ) 
dx 
a) 
Multiply imerator and denominator by proper factor } 
—)2 
“— NIix+X 
/ at PSD 
, ix 
7-2 where X= X1+X9...X, 
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= Zy 
- | (y-X)e dy 
toe - X 
“en 2 
“Sy 
€ dy 


X1+Xor : Xn 





n 


Thus in the case of a Gaussian distribution, the function (g) reduces to the average 


and no improvement over the usual practice is foynd. 


If the distribution is rectangular, then 


n 

p=A S(X,) S(X%9) S(X,) 

where S(X) = 1 |x| <a 
and S(X) = 0 |x| > @ 


and A is the normalizing factor 


Therefore, from equation (4) 


5 he S(X+X 1) S(X+X 9) S(X+Xn) dx 


gs 
F im S(X+X1) S(X+Xo)...S(X+X,) dx 





Let X; = the variate having the lowest value and Xp 
value 


Because of equation (7), the integrands vanish for 
X+X;< -a and X+Xp,> a 
or for X< -a-X, and X <a-Xp 


Hence, equation (9) becomes 





a-Xh 
- X dx 
g- -a-X, 
a-Xp, or 
dx 
-a-Xy 
Xhn+X 
g hed 
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(7) 


(8) 


the variate having the highest 


(9) 


(10) 








Equation (10) shows that for a rectangular distribution, the average of the highest 
and lowest variate is the improved method of calculation. 


For some distributions, equation (4) yields expressions considerably more compli- 


cated Hence, for Cauchy's distribution 


\e 
P  . 2 2 
(1+X 1“) (14K)... (1+Xq*) 





gives for 


8(X 1 +Xo+X3) -~6X 1X 9X3+X 1 7(X2+X3) +X27(X1+X3) +X37(X1+X9) 





. 2 2 2 
2(Xj +X9 +X3~ -X Xo -X,X3 -X9X3 +12) 
APPENDIX 
Given g(X1, X2,...Xn) = g(KitU, KX2at+U,...Xnt+U) -U a) 


and t= / g?pdxi dxz (2) 


where p = p (X1,X2,...Xn) 


which satisfies equation (1) and which minimizes the right 


To find the function (g) 
(1) and (2) to the new set of variables given 


side of equation (2) transform equations 


as follows: 
xX xX} 
Y} X9-X} (3) 
y2 = X3-X) 


F(X, y1, y2-.-Yn-1) = F(X+U, y1, y2...yn-1) -U (4) 

where F(X, yl, y2...yn-1) = g(%1, X2...Xn) (5) 
Equation (2) becomes 
.2 

I - | F? g dx ay) dys... dyn-] (6) 


where §(X, yl, y2,...¥n-1) = p(X1, X2,...Xn) (7) 
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Since equation (4) must hold for all values of U, if we expand the right member of 
equation (4) into a Taylor series in U about the point X,yj,y2...yn-1, we obtain 


oF _, (8) 
ox 
or F = X+f(yl, y2,...yn-1) (9) 


where (f) is any arbitrary function of the y’s. Since equation (9) was a necessary 
condition on the solution of equation (4) and substitution into equation (4) shows it to 
be a sufficient condition, equation (9) is the general solution of equation (1). 
For (I) given by equation (6) to be a minimum, it is necessary that 51 = 0 
or dI = 2 {¥F¢ 5 Fdxdy; dy2.. .dyn-1 = 0 (10) 
for all dF compatible with equation (9). 
But by equation (9), the form of F is determined except for the form of (f), hence 
dF = df (11) 
Substituting equations (11) and (9) into equation (10) yields 
dI = 2 f x+1) $ Stax dy; dyg...dyn-1 = 0 (12) 
Since df is a function not dependent on X, equation (12) may be written 


2 [bt ayy dye... .dyn-1 /Ox+t) pdx 2s (13) 


Equation (13) may be satisfied for arbitrary (df) only if 


[ext g ax = 0 


or f -} X9 dx (14) 
dx 
F x] -J Xp ax a aes [ (x-xX}) dx (15) 
Jt& dx 


If we transform the right side of equation (15) back into the old variables, we ob- 
tain 


- 
— J xpeeexs, X+X2,...X+XKn) dx 


s* = (16) 
p(X+X 1, X+Xo,...X+Xn) dx 





Since (g) given by equation (16) is the function of minimum variance, the function 
(g+a) must have greater variance than (g) for all values of the parameter (a) different 
from 0. This can be true only if 


g = 0 (17) 


Hence, it follows from equations (1) and (17) that if the distribution (p) is given ex- 
plicitly with a location parameter 6 
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p = p(X,-8, X2-6,...X 


n~9) 


Then (g) given by equation (16) is an unbiased estimate of the parameter 6. 
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The Differential Equations of Structures 


Clayton J. Pajot . 


Int roduction 


The following discussion is intended to show the derivation and application of some 
of the differential equations most commonly used in structural engineering. 


The discussion is subdivided into the following headings: 


1. Derivation of the fundamental equations relating to bending moment, stress, de- 
flections and the beam constants 


2. Application of the fundamental equations to beams, columns and beam columns. 


3. Derivation and discussion of the energy method of finding deflections (castigliano 
theorem) 


4. Extension of the fundamental equations to certain cases of plates or slabs. 
Assumptions 


The assumptions made here are those usually made when discussing problems in 
elementary strength of materials, namely that: 


1. The material is perfectly elastic (stress is proportional to strain). 
2. The material is isotropic. 
3. The principie of superposition holds. 


4. The ratio of stress to strain is the same in tension and compression. 


REFERENCES 





Timoshenko - Strength of Materials - Vols. I and II 
Timoshenko - Theory of Elastic Stability 
Niles and Newell - Airplane Struct ires - Vol, Il 


PART ONE 


Pp 
2 48 











a — = 
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FIGURE 18 a b 
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Let Figure 1A represent a beam of length (L) and loaded as shown with concentrated 
loads Pj, Pg and P3, and a uniform load of (q) pounds per unit length. 





1 
Orientation Lecture delivered before the Society, March 31, 1950. Clayton J. Pajot 
is Chief Stress Engineer, Nash Research Department, Nash-Kelvinator Corporation, 
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Figure 1R represents the deflected position of the beam in exaggerated form. Evidently 
the upper fibres of the beam are compressed or shortened, whereas the lower fibres 
are elongated. There is then some plane between the upper and lower surfaces of the 
beam in which the fibres are neither shortened nor elongated. This is called the 
neutral plane of the beam and is represented by the broken line passing through the 
points (a) and (b) in Figure 1B 

















FIGURE 2A FIGURE 28 


Figure 2A represents an enlarged view of the element a - b. The external moment 
is shown to the left of the figure and represents the moment obtained from the ordinary 
moment diagram for the beam. It is assumed to be positive when it causes the beam 
to be concave upwards as shown. The triangles on the face (b) represent the variation 
in stress onthe internal fibres of the beam. The stress is assumed to be proportional 
to the deformation and therefore has a zero value at the neutral plane a- b, a maxi- 
mum compressive value ao, at the top and a maximum tensile value o 9 at the bottom. 
The stress at any distance (c) from the neutral axis is represented by o. 


Figure 2B shows a section of the beam at the point (b). For this discussion it could 
be of any shape symmetrical to the vertical axis. For ease in drawing, it is assumed 
as shown. The line b - b is the intersection of the neutral plane with the vertical 
section and is called the neutral axis of the section. 


(dF) is the total force due to fibre stress acting on the differential area dA. 


THUS OF =o: dA 
FROM SIMILAR TRIANGLES O.C OR o= 2 ¢ 
% Ce Ce 
we dF = Ce cd , 
Ce 
AND Co Ce (1) 
F = dF = ce € dq = Ce cA 
G e Je, Ce 


Now (F) represents the total force acting on the face (b) of the element. Since the 
element is in equilibrium due to moment (M) on the face (a) and the forces on the 
face (b), and since (M) represents a couple which has no resultant force it follows 


that Ce 
[o- O oR & CA =0 
¢, 


Ce 
gut oO AO and Co#°O 


therefore @ = o and b - b is the centroidal axis of the section. 
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Let (dM) be the moment due to the differential force (dF). Then the external moment 


(M) is equal to the total moment due to the differential forces along the face ¢b). 


Ce Ce 5 . 
OR M = dM = cdF - 2] cdg 

Ce c 
c la 2 
' P ¢, 2 

I= oda 
¢, 
where (1) is the moment of inertia of the area of the section about axis b b. 


Oz om 
THEN M= Ce . 2 ef 
MC (2) 


ND r= — 
I 


In Figure 3 i hown the element a - b of the beam in its deflected position. (o!) is 
the center of curvature and 7 the radius of curvature of the plane a b. 
The line (bf) is drawn parallel to (ad) Then (df) is very nearly equal to (ab) and (fg) 


is the elongation of an element of original length (ab) 


Let € = unit strain at distance (C) from the neutral axi 


4E te 
dt ab 





FIGURE 3 


then Es 


considering the similar areas (bfg) and (abo!) 


J 
t6 


c 
abl ha f 





Let (E) the modulus of elasticity of the material \ 
-_ sTRESS_ o -— / ) 
C = STRAIN = = OR aie in Be 
7 \! 
Cc Ec -fa/ / A 
THEN — = — OR c= — at 
? -° “of tA g 
ee - 
Me —~ i 


ae: a. ar ae (3) 
fp T pP EI 


{ \T+4T 
ca ae Te 


FIGURE 4 


Figure 4 represents any curve Y = f(x) whose curvature 2’ \/P isto be determined, 
From the figure As 3-At (approximately) and 1/7 At/As (approximately). 














; Se 
THEREFORE " BP ~ a8eo 4s ds 
-/ 
IN FIGURE 4 oY =tan T OR ro = 
dy dx 
d@y / 
| I , ALSO dS. | 7 + (ex \* 
_ —1d¥ dx 
dx/ dT d*y 
dr, dk _ __ dee = < 
AND ds * ds [ 2 (gr) ;? r 
dx d 
M SF 
x2 
—_———— = = "a . 
"EI f,, (x)? |* - 
L \ dx : 


Since, in the case of beams, the slope of the neutral plane at any point is in general 
small, (dY/dx)? in (4) may be disregarded when the formula is applied to members 
having small deflections 


mM_ d*7 
Then (4) EI = Jdxe (5) 





which might be called the fundamental differential equation of structures 


PART TWO 


(A) Beams 
Let Figure 5 represent a beam having a uniform load of (q) pounds per unit length and 
the concentrated loads P); and Pg and having end reactions of R) and Ro The dif- 
ferential equations of deflection will be written for the beam. 




















O @ @ 
—_— ssi a —=————— - 
i. cnad t 
1GUF 
R FIGURE 5 he 


Note that the concentrated loads constitute points of discontinuity for the moment 

equations and that therefore three separate equations must be written, one for interval 
(> , 

i), one for interval (2) and one for interval G) 








The moment equations for the three intervals are: 











Mm = RK - $€ o<xZa 
Me = RX - XP _ PCx-a) a<x<b 
Ms = RX — -- P(k-a)- Fe (x- b) b<x <i 
Using equation (5), the corresponding differential equations are: 
ty rv, iT R x= 
a - / = x - 6 
dx ~ Er ~ etl ' Z | o<x<a (6) 
ey , 2 
fae = Me — | &x - 4x - P(x-a)) a<x<b (7) 
ox Ez EI z 
Jfy . Ms ; | $X* _ prx 
= =——| &Y¥ - - P(K-a)-Blx-b))  b<xKeL 
dx@~ ef EFI : é (8) 
Note that if the cross section of the beam is not uniform throughout that the moment 
of inertia (I) will not be constant This condition may be expressed as 
I = Fix) 


(= Comes 


Let Figure 6A represent a column fixed at its lower end and with an eccentric load 
(P) at a distance (e) from the center line of the column. Figure 6B represents the 
deflection curve of the neutral axis of the column and Figure 6C shows the deflection 
curve rotated to conform to the coordinate system as previously chosen. 


pP s 424 Y 
tee . 


4 
FIGURE 6C 














Yo 











Mp! skh 
FIGURE 6A FIGURE 68 


From Figure 6C the moment equation is: 


M = Pie +d - Y) 


Then, from equation (5) 


Ps ? o 
f@TiM.L(e+d-y) = k*(e+d5-y) 
dx= Ez EI 

2 
WHERE Kael 


ELT 


OR dy + Key = K*e+d) (9) 
dx® 
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The solution for (9) is 


gf = e( sec Ly =f) (10) 


Equation (10) is plotted in Figure 7 


A ASYMTOTE B 









$a e (sec Lyf -1) 








FIGURE 7 
If there is no eccentricity, e—-o. Then the so-called critical or buckling load is 
P., where P., = lim P. 
e@ —e 


The value of P.,. is given by 


Pp 


ce 


z 

= TEE (11) 
Fl 

This is called Euler’s column formula and corresponds to the condition shown in 

Figure 8 


Theoretically, when a concentric load (P) is placed upon the column and then gradually 
increased, the column should deflect in accordance with the curve oAB of Figure 7 
which means that when the critical value of (P) is reached the column should suddenly 
fail. In actual test this failure does not occur, but when P., is reached, it deflects 
somewhat as shown in Figure 8 


FIGURE & 


The reason for this discrepancy can be attributed to the fact that when this condition 
occurs, we are dealing with large deflections but equation (6), upon which our solution 
is based, is to be used for small deflections only 
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For large deflections, equation (4) should be used. Applying this equation to Figure 


6C we have 


ga? 


dxe - Mw -P « ¥ 
[rave | f= £ (e+s-¥) (12) 


i} (= 
i. K / 7 








In the case of a column with hinged ends, the solution of this equation has been found 


to be 
§= 446 ze -1|1-£(E -1) 8 


Figure 9 shows the graph of equation (13). It will be noted from this graph that when 
loads above the critical are reached that large deflections occur with small increase 


in load 





+ 


j 
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FIGURE 9 
BASED UPON LARGE DEFLECTIONS FIGURE 10 


Figure 10 shows the manner in which a column having hinged or pinned ends deflects. 
This type is usually called the fundamental case The critical load for this type, 


based upon the equation for small deflections, 1s 


P 2 
pp. war (14) 
CR ——— 
sé 
<a 
A condition which sometimes occurs in practice is that shown in Figure 11A in which 


an initial eccentricity (e) exists before load is applied. The column deflects into the 


form 11B with a deflection § as load is applied. 


P 


‘a 


lad 





FIGURE 11A FIGURE 11B FIGURE 12 


Another type of column sometimes encountered in mechanical and aeronautical engineer- 
ing work is the stepped column shown in Figure 12. 


~] 
~] 








An interesting discussion of both types may be found 
issues of Product Engineering magazine for the year 
Meier 


“Uniform and Stepped Columns” and is by J. H 


in the October and November 
1949. The article is entitled 


Note that in all this discussion of columns it is assumed that we are dealing with long 


columns; that is 
pared to the length of the column 
columns 


(C) Beam Columns 


Beam columns occur frequently in aeronautical work. 


columns having cross-sectional dimensions that are small as com- 
These are sometimes referred to as slender 


A beam column is simply a 


member having some form of transverse load in addition to its axial load. 


Figure 13.shows a beam column having only one 


into two intervals or parts numbered as (1) and (: 


transverse load. It is subdivided 
on each side of the transverse 


load, as the load constitutes a point of discontinuity for the moment equations. Two 


equations are then developed, one for each part 
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FIGURE 13 
M,=- - QBAX _ PY , M, = @2(4-x) _ P% 
L Z 
FOR FIRST INTERVAL O< X LBW 
J*y¥_ mM __ QbX _ py 
dx® ET EIL EI 
Lf «a 
LETTING ET 
oy. 0 « «~ is 
dx® EIL (15) 
FOR SECOND INTERVAL a<x<l 
d*y_ Me_ _ Qa(tL-x) PY 
dx* EI EIL EI 





2 
AND LETTING -fL = K 


&£I 
2 

d7Y 4 K*Y = _ _ Qa (.-x) 

dxé@ EIL (16) 
THE GENERAL SOLUTION OF (15) IS 

Y= Asn KK + B cos Ky - QbxX (17) 
PL 

Where A sim Kx + Bcos Kx is the complementary function and - a is a particular 
solution PL 


Similarly, the general solution of (16) is 














Y% = C sn Ax + D cos Kx —- Qa (l-x) (18) 
SINCE y, =O WHEN xXz=O PL 
%=0 WHEN X=Z 
AND = Ye AND GY d¥% WHEN X=Q@ 
ax dx 
The constants A, B, C and D may be evaluated 
A= Q@ sm Kb , Bz=0o , Cz_ Q sin Ka 
PK sim KL PRK tx AL 
AND D= @ sin Ka 
PK 
Then, substituting these values in (17) and (18) 
Y= Q smkb smkx ~ Q6xX o<xXéa (19) 
PA sin AL PL 
a<x<b 
=-~Qsnka sinkx 4 @ sin Ka cos kx . Ya(L-x) (20) 
PRK tan AL PL 
Substituting the values of Y; and Y2Q in the moment equations: 
M, = — Qbx _@Q sunkb sinkx és Ge x O< xXx <a&@ (21) 
4 K simKLl 
Qa (i-x) 
a“ 
Mp= @a(L-x), Qsnka smkx_ car Ka cosKx + (22) 
é K tan KL bexel 
Ymax may be determined from (19) and Mmax from (21) or (22), This method may 
be extended to a condition having several concentrated loads. Note that 
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Another type of beam column is shown in Figure 14. In this case, in addition to the 
axial load, a uniforn distributed load of (q) pounds per unit length is applied, a 
moment or couple M, at the left end and a moment Mo at the right end 
& 
” “2 
== 
M, soatintitiam oth } Mp 
PTT Th 7 Wh 
Litt tii} 
Lime 3 amen aa oY length 
Fe 
P - ee, fe 
M, ‘¢ Pn Jt = + a —_ ”) Me 
_s £ — 
FIGURE 14 
The method I n tt iS€ ight different than in the previous case 
Although a differential equatior he same form is derived 
0 = #4 1, - We £L=- 24 My - Hh, 
— # en eS —— + 
” L— 2 £ 
f D g "fs % 
2 
2 a -£- Pf SZ 4 gu 
/ x ,, @% 
2 / 
g 
GH... - y + @ 
ox? a x* 
2 
BUT FROM (5) a ae 
dxé £Z 
2 
AND THEREFORE Le. te. Z 
p " ox? £L 
LET —_—_— 
iil 
THEN /*A7 2 
aX 
THE GENERAL SOLUTION | 
M2 AsnaKXK & B cos KX *# wr) (25) 
where A sin Kx + B cos Kx the complementary function and q/K? is a particular 
lution 
WHERE 4 Sin ek + £ Cos A \S THE COMPLEMENTARY FUNCTION AND 
BYe2 IS A RTICULAR SOLUTION. 


WHEN X =O, Af m= AZ AND WHEN vu Z, Af 
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a 4 
aa hn a s Zaiteoe & 
AND ‘7, ge 





° ~/77 x & 7 
Substituting the values for and in (25) 
4 
Yo — 42 ~-(M4 - F/ C05 KZ 
M= ee Gwen Spe Kx 4 (i - =, coskxe £ (26) 
If there are no end moments Mj= Mg? 0, then 
Z e 
- > + COS a 
2 2 
M=x —<~ £ . Sew. newtr «2 
Sin KL , 4 yal 
g - 
¢ a & =-/ < 
_ # 2 Set) Sm KK - COSKK ¢ | (27) 
ge L S10 Py A 
Mmdx may be obtained from (26) or (27) The equation of the deflection curve may 


be obtained by substituting equation (26) in (23) and then solving for Y9. This gives 





- ~ _/ P 
Yehitt-Le, E* Me — qa ~ (hs ~ ga) 008 KL 
DP | p ‘ a olen —- Sia KX 
Sioa KL 
_(M,-4% \cos kx 2] (28) 
Ke, Ke | 
The slope equation is obtained by differentiating (28) It is 





g. ) 
oe afk _ ek ygx — ene FoJeos kl 
[ Sia KZ 


KoosKx 


Beam columns continuous over several supports may be solved by applying equations 
(26) and (29) to each span independently 


A value for either the slope or the moment must be known at the beginning of the first 
span and at the end of the last span It must be remembered that the slope at the 
end of any given span is equal tothe slope at the beginning of the adjacent span. Thus 
a series of simultaneous linear equations may be set up by means of which all the 
moments over supports can be obtained Then, knowing the end moments for any 
span, the maximum moment for the span may be obtained from (26). 


PART THREE 
(A) Energy in Beam Based Upon Bending Moment 
Let Figure 15 represent a beam of length (L) with concentrated loads and uniform 


loads as indicated. The element a - b is any element of differential length (dx) located 
at a distance (x) from the left end 
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! -. - 
[ 445 per wat length 
| wns 
a x _.!|_2* 
| 
_——_—— - é 
R, FIGURE 15 ‘2, 


Let (M) represent the moment at point (x) onthe beam, obtained from the usual beam 


moment diagran 


Figure 16 represents an enlarged view of the element a - b shown with the deflection i 


due to the moment (M) exaggerated 
dU work stored in element a - b due to bending 
=. Add ‘ 
2 * O 
WHERE og a ANGLE OF BEND 
FROM FIGURE 16, 


SX =a lAb OR <a. 
,-° - 





USING EQUATION (3) 


THEN THE TOTAL ENERGY OR WORK STORED 
IN THE BEAM IS GIVEN BY 





, 4 
X —-/cSH= (4X (30) 
FIGURE 16 


Then the total energy of work stored in the beam is given by 


(B) Castigliano Theorem 


Consider a beam loaded as shown in Figure 17 
p | |6 Rk |f 
= | ! . = 








| gf 4% per vai? Jeng (bh 
SSS 


x 


zB, 
FIGURE 17 & 








Let (U) be the energy stored in the beam due to the application of the loads shown, 


Then U=f(q, Pj, Po, P3 ----- P,----P,). 
Suppose that at point K an additional differential load dP, is placed (See Figure 18). 
2 
Pp 2 B R RB 





| g (bs. per PT 4eg 7b 








4 


2, FIGURE 18 R 
2 
Let (U') be the total energy now stored in the beam. 


+ a 
Ys wa tk 





Now consider the case in which the differential load dP, is first placed onthe beam, 
and then the given load as shown in Figure 17 is superimposed. The total energy now 
stored in the beam is the same as in the first case and may be again represented by 
(U'). Thinking of (U') from the standpoint of this second procedure it may be repre- 
sented as follows: 


where (U) and dP; are defined as before, but ¢, is the deflection at (K) due to the 
application of the loading condition represented by Figure 17. 





Since the two values of (U') are equal 
“u+ EIR =K+d2) ds, 
oY r 0X 
OR ay A JR md AR OR o = ok (31) 


which is the theorem of castigliano 


From equations (30 and (31) an expression for the deflection may be written as follows: 


6-%.2-( SH. PM Jx (32) 
- ae ASR, 2ZETL LL £2 Akh 


An expression similar to (32) may be written for curved members or beams, the re- 
quirement inits application being that tne radius of curvature of the member is large 
as compared to the depth of the member. Figure 19 shows a curved beam or arch 
and equation (33) is the corresponding equation. 








BR 
a /® £ = Length of 
Me ——==> Nevire! Axis 
P Cox \ , 
y \ 
FIGURE 19 J 
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f WM oo M Js 


Y, ~ J fr I ~ (33) 


} 


In equation (33) (ds) is a differential distance measured along the neutral axis and 
S;, is the deflection of the member at point (K) and in the direction P, due to the given 
aading ndition 


(C) Discussion 


Equations (32) or (33) may be solved mathematically when (M) is expressed in terms 


the known loading and (x) or (s), provided that the moment of inertia (I) is either 


4 constant or can be expressed as a function of (x) or (s) 


However, either equation (32) or (33) lends itself very well to graphical analysis when 
(I) cannot be expressed satisfactorily as a function of (x) or (s) or when the loading 
condition is such as to make it difficult to express the moment mathematically 


Although most of the equations that have been discussed here are related to deflections, 
deflections in themselves are not necessarily the most improtant thing to consider 
when dealing with engineering problems In many cases the object of obtaining de- 
flection equations is to use them in finding reactions, and ultimately moments and 
stresses in the indeterminate members or structures 

For example in Figure 20A is shown a curved member represented by the curved 
line shown with both ends (A) and (B) fixed, and with the loads indicated 


Figure 20B shows the same structure but the fixed support at (A) has been replaced 
by the three reactions Ha, V and Ma. In this figure the loads Pj, P95, and P3 are 


2 


known, whereas the reactions are the unknowns to be solved for 


P 
TZ i"2 
a t 


{7 IIA 


Aa Mg 
z, ) FIGURE 208 777 


“4 
Let M f (H, Va, M,, loads) where (M) is the moment at any point such as (C) 
ind is obtained by taking moment f all forces to the left of point (C) 
} eq it I { 


dij =O (34A) 
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where §} is the deflection of the end (A) in the direction H, and is zero due to the 
fact that end (A) is fixed 


Similarly, 





z 
- M oM 
O J et Vp é3 “0 (348) 


where § y is the deflection of end ) in the direction V, and 


"— @i so Ac 
a Fr on fo) (34¢) 


where 6, is the rotation of the end (A) 


Equation (34A) will be of the type 


iy, 4%, 4,, leads) =o (35a) 


4 
and similarly (35B) and (35C) will be of the form 
4, Cy Vy, ty, leads )= o (358) 


j3 4 (Ha, V, , My, joads ) = Oo (35¢) 


Equations (35) are simultaneous linear equations that may be solved for the unknowns 


H,; Va, and M, 


The reactions at (B) may then be determined from the equations of statics and also 
the maximum moment and maximum stress in the structure 
PART FOUR 


(A) Rectangular Plates 


Let Figure 21 represent a thin rectangular plate of thickness (h). M, is the moment 


per unit length applied along the (Y) axis as shown and M,, is the moment per unit 
length applied along the (X) axis : 











; aay 
e jun” 
il 2 a 2 
[ a f { 
| M, 


FIGURE 21 


(abcd) is an element of the plate of dimensions dx and dy shown enlarged in Figure 22 
Compare with Figures 2A and 2B 


Let o, be the stress per unit area on this element in the (x) direction and vy the 
stress per unit area in the (y) direction 
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FIGURE 22 
(a'b'c'd') is an element of Figure 22 of thickness d z and is shown separately in 
Figure 23 with curvatures due to the bending moment The points (O') and(O") are 
the centers of curvature in the planes parallel to the x 2 and y 2 planes respectively. 
z2 
~=.= => AND é, - 
F i‘? 


where €& and €éy are the strains in the (x) and (y) directions respectively. Compare 


x 
with derivation of e ®/po from Figure 3 ' 


.O 
/ 
A — —=— — i Sr" 
ate ~ EG A Oy } 


/ « 
AND < — (o_- Oo 
x 
where is Poisson's ra 
a 4 
Eliminating € , from the last two sets of equations; d 4 
£2 4 / 
Ge Bear / a tA 5 
al (2 (jy 











FROM FIGURES 21 AND 22. FIGURE 23 

rf , 
M dy -/ 2 Wde oR Mw “7 ea de 

-# oe 
THEN MW j? e* ( 7 + M , d2 
= re 2 — 
‘ 44 — (x (y 
Zz 
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LET > =z (36) 


a / 
THEN A == / Md 





('y 
AND, SIMILARLY, Ad =e D/ 4 + he = (38) 
[ 2 2 
AS FOR EQUATION (5) 
| ee ou“ = “ar 
a ~ ~ ox ~~ -—_— wa 


where (w) is the deflection of the plate and the negative sign is introduced since the 
Z axis is taken as positive downward 








3 " O*ur O2ur | 
THEN Ad 2D, see + i ore) (39) 
iad: i Chis fur (40) 
_ oe of 5 yé, a ox- 


Compare equations (39) and (40) with equation (5) 


(B) Circular Plates 
In Figure 24A is shown a thin circular plate of thickness (h). The loading on the 
plate is symmetrical to the center and is the same in any section through the center 
as shown in Figure 24C (R) represents the reaction along the circumference of the 
plate and is in pounds per unit length. Figure 24B shows the curvature in any diametral 
section and also the radius of curvature /, inthis section. Figure 24E is a diametral 
section but shows the radius of curvature 75 ina direction perpendicular to the section 
(circumferential direction). Note that the center of curvature corresponding to @9is 
on the axis O Zas elements such as the one through (A) form a conical surface with 
the apex at (O") The moments M, and Mg shown in Figures 24B and 24D are mo- 


ments per unit length of direction ab or bc. 

dur 

From Figure 24B —— 
e Pp a x 


Note that (w) is taken as positive downward hence the negative sign 





2 
sleee ton ae a= is SEE EQUATIONS (3) AND (5) 
- AX gx 
FROM FIGURE 24€: = Jv w= 
G fA ingmp 
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Jo | 
\ 2 crt FIGURE 24D 
+ 
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Mz 
\ 
“i lA 
FIGURE 24A 
a 
\ 
wie 
Lm idg 
| p — xX 
M, 
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ewan os yr “ 
A? o’ | | 
| |_* atx | 
2 
e / 
O4 =f? 
a“ 
“ee & I TX eo teat 1 
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oO’ 
pe RADIUS OF CURVATURE 
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DIRECTION 
FIGURE 24¢ 


88 





Pla S$ af 47 — Mz 














sie x a 
WA FIGURE 25C 
SM, 
~+—— JX = 
| va’ l - ins ad=xde 
O} \SLZA) Tt ——* 2 =tenende 
. @ 4 Ja 
| Dp a SX 
| FIGURE 258 ox 
2 
Substituting in equations (37) and (38) 
MZ +1KG) = , ¥) (41) 
(7A 2) =4 Ger 
/ / J 
M,= D(> fr rire LP . 
. (ZB AP x ~ ox $8? 
Taking the summation of moments about point (A) in Figure 25B 


2Mg=(M,* Te Uk) (x4 dx ee. seen -2M, a2 Sx 
+ Qr @ dn) (« ¢ bx} d0dx = O 


dividing by d@ and expanding 


M,x + on £4 dx + M,dx t ae, -M,x - M, —— 
ge ax)? 4 Qlax)* a So (ax) “= O 


Combining terms, dividing by dx, and dropping differentials of higher order 


M,+x L4-4,7 Qx =0 (43) 


and substituting the values of M, and Mo from (41) and (42) into (43) 


i 


a 
(gh + ut) +x0f Fh re 3 -) 


- (Fras) 4+ @X=0 
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2 
é 7“? Gf vwh.j.§§Lastha= 4. 
OR ax tA TX Gat ax A x x “A ax ra) 
Jd + —— Q | 
—_<——— <= —_—_——_- — 44 
AND me? X - x2 D (44) 
WHERE PD =  —_ SEE EQUATION (36) 
Ji C0 /-A4) 
For a circular plate with a load (P) in the center 
P=Z2TxQ Rn G= ” AND (44) BECOMES 
277X 
2 P / 
Chest -GBhea «~~: (45) 
dx* K dK x 27D x 
For a uniform load of (q) pounds per unit area 


7X*g =27XQ OR Qu Xx AND (44) BECOMES 
i te. theidis 


Jx*# X dx x? 2D (46) 


Note for integration purposes that 


f6,1d@ @od [st a ] 
ax” X dx x2 dx lx gx (#Y| 
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THE THERMAL DESIGN OF EVAPORATIVE CONDENSERS 
Hugh J. Scullen! 
Abstract 


As a practical example of a method of analyzing and using heat transfer tests the 
data from an experimental evaporative condenser are first reduced to individual film 
coefficients Then these are inserted in general equations which can be used to 
estimate the thermal performance of any design. 


The purpose of this paper is to show how fairly elaborate mathematical methods of 
analysis and synthesis may very profitably be applied to test data of a relatively low 
degree of accuracy when the object is a comprehensive rather than precise picture 
of the character of a heat transfer device 


The particular kind of device used as an example is what is known as an “Evapor- 
ative Condenser.” This is essentially a forced-draft, air-cooled condenser provided 
with means for maintaining a layer of water on the air side of its heat transfer sur- 
face so that heat will be removed from it by evaporation as well as by convection. 


The heat dissipated by an ordinary dry-surface, air-cooled condenser can be com- 
puted by the well known equation 


q U ot 
At] eq. |! 
I 


where the factor U, called the “Overall Coefficient of Heat Transfer” is a function 
of the air velocity and, to a small degree, the temperature and rate of condensation. 


It is common practice to show the results of air-cooled condenser tests by plottin 
if 5 


this coefficient U against the air velocity V, neglecting the effect of temperature 


and rate of condensation A much more useful way of summarizing test data is to 
give equations or charts for the values of the individual film coefficients which com- 
bine to form the overall coefficient The relation between the coefficients is 

l l l 

, Ct SS > 2 

l he Rh, eq. « 


As shown in detail in Appendix 1 of this paper, the heat dissipated by a wetted-surface, 
air-cooled condenser can be computed in an analogous way. 





This paper was presented in abridged form at the Semi-Annual A.S.M.E. meeting 
held in Detroit June 17, 1946 and at the first meeting of the Industrial Mathematics 
Society, held February 28, 1949 The data and some of the figures are used with 
the permission of the Nash-Kelvinator Corporation, for whom the tests were run in 
1935. Hugh J. Scullen is a consulting engineer in Detroit 








REGULAR MEETINGS 





The first public meeting for the purpose of organizing a society to promote the use ' 
of mathematics in industry was held on February 28, 1949. Robert Schilling, H. J. 
Scullen, J. E. Freer and A. W. Jacobson participated in the discussion on “Industrial 
Mathematics” with E. | tiopelle as chairman The meeting instructed the or- 
ganizing committee to proceed with the preparations to formally establish the Industrial 
Mathematics Society The constitution was adopted at a meeting on May 24, 1949, 
and the officers were elected Dr. R. S. Burington of the Bureau of Ordnance 
spoke on “The Importance of Anchorage in the Application of Mathematics to Industry.” 
After the preliminary meetings, the following four regular meetings of the Society 
were held during the fiscal year of 1949-1950: 

October 14, 1949. Speaker: Hillel Poritsky, General Electric 

Company ~ Sut ject: “Some Industrial Applications of Mathematics 

December 2, 1949 Speaker: A.D. Moore, University of Mich- 

igan ~ Subjec t: “Fluid Flow Mappers 

February 14, 1950. Speaker: A. W. Jacobson, Wayne University 

Subject: “ Analogue Computation and Application to Engineering 

Problems.” This meeting was followed by a two-day conference 

on the subject in which E. A. Caldwell of the Massachusetts 

Institute of Technology and a large number of local engineers took 

part 

May 24, 195¢ Speaker: Irvin Travis, Burroughs Adding Machine 

Company ~ Subjec t: “Digital Computation - Its Business and En- 

gineering Applications Two conferences were held the next day, 

one on “Logic, Use and Programming” with A. W. Burks of the 

University of Michigan, and the other on “Design and Construction 

of Digital Computing Machines” with C C. Chambers of the 

University of Pennsylvania 

ORIENTATION LECTURES 

Three orientation lecture arranged by the Group Activities Committee, were held 
during the year These lectures by Chenea, 
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Hay and Pajot are being published in 
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STUDY GROUPS 
The following five Study Groups were active during the year: 
Gear Design: 


Group Leader: Douglas Hughson, Nash-Kelvinator 
Adviser: W. F. Vogel, Wayne University 


Computing Machines: 


Group Leader: R. A. Roggenbuck, Ford Motor Company 
Adviser: A. W. Jacobson, Wayne University 


Servomechanisms: 


Group Leader: Jorma Sarto, Chrysler Corporation 
Adviser: W. H. Bixby, Wayne University 


Electron Circuit Theory: 


Group Leader: P. J. Willson, General Motors Corporation 
Adviser: K. R. Symon, Wayne University 


Stress Analysis: 


Group Leader: Egon Benesi, General Motors Corporation, 
who also acted as an adviser. 
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SUGGESTIONS TO AUTHORS OF PAPERS 





An author wishing to contribute a paper to “Industrial Mathematics” should submit 
an original and copy of his manuscript, typewritten, double-spaced on 8-1/2 x 11 
in. paper Compliance with the following suggestions will facilitate the publishing 


of the paper 


HEADING 
] Title of paper 
2 Author or authors’ names with academic degrees 
3 Present affiliation of author or authors 
(a) Position held in company (or university) 
(b) Name and address of company 


B. ABSTRACT 
Include a short summary to be printed with the paper 


REFERENCES 
l Should appear as ; l at e end of the paper 
2 Form 
the following order: authors’ last names 
le, name of journal in which it appears, 
volume ni “r, page number and year 
Example: 1) Jones, J., “The Fundamentals of 
Vol. 153, p. 592 (1936). 


book -give in the following order: authors’ last names and 


Physics,” Jour- 


title of yk, publisher, city in which it was published, 


ear of publication, and page reference 


Example: (2) Janes, M. E. and Mason, E. R., “The History of 


Nuclear Physics, The Atomic Book Co., Inc., New York, (1948) 


D ABBREVIATIONS 
Follow the America: entative Standard Abbreviations for Scientific and Engi- 
American Standards Assn. 


FIGURES AND FIGI 


t larger than 8 x 10 in 
numbers only should be written lightly on the back of the 
a = 
pencil Captions should be typed in order on a separate 
r and should contain the figure number. 
Drawin 
(a) Mi be in India ink, preferably on tracing linen or tracing paper, 
eptable on good quality heavy white paper. 


re not 


o be included as part of the drawing but treated 
photographs 


The editors reserve the right to make minor editorial changes in the manuscripts, 
so all paper lay conform to the over-all rle of the Journal 


rhe author will be supplied with 20 unbound copies of his paper, and additional copies 
may be obtained at the author’s expense These should be ordered when the manu- 
script is submitted to avoid the necessity for a rerun 





